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PREFACE 


This  Report  is  a  product  of  an  investigation  of  border  infiltration 
undertaken  by  The  Rand  Corporation  for  the  Advanced  Research  Projects 
Agency*  It  completely  discusses  the  logic  underlying  PREMIR  (a  Predic¬ 
tion  Model  for  Infiltration  Routes)  and  provides  information  on  the 
model  and  illustrative  results  from  its  first  application.  Most  studies 
of  infiltration  deal  with  the  problem  in  the  context  of  detecting  enemy 
logistics  and  operations  through  aerial  reconnaissance,  prisoner  or  in¬ 
former  interrogation,  reports  from  the  field,  etc.  These  studies  have 
generally  not  been  concerned  with  explicitly  developing  the  relationships 
between  patterns  of  infiltration  routes  and  the  infiltrators1  criteria 
for  choosing  them.  In  sharp  contrast,  the  present  inquiry  is  intended 
to  more  fully  identify  and  formulate  the  essentials  under  which  infiltra¬ 
tion  routes  may  be  systematically  analyzed,  modeled,  and  predicted. 

Although  the  analytical  techniques  employed  in  this  study  are  well 
established,  the  particular  application  of  them  in  an  investigation  of 
border  infiltration  is  believed  novel.  For  that  reason,  the  information 
presented  here  should  be  of  prime  concern  to  others  interested  in  the 
problem  of  infiltration.  But  as  the  study  illuminates  certain  aspects 
of  infiltration  that  have  so  far  received  little  attention  in  the  liter¬ 
ature,  it  should  also  be  of  interest  to  a  much  wider  audience. 
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SUMMARY 


This  Report  investigates  some  conditions  under  which  infiltration 
routes  can  be  systematically  modeled,  analyzed,  and  predicted,  and  iden¬ 
tifies  some  areas  where  the  resulting  information  may  be  efficiently  uti 
lized.  PREMIR  (Prediction  Model  for  Infiltration  Routes)  is  primarily 
a  mathematical  computer  model  which  by  means  of  established  analytical 
techniques  arrives  at  a  set  of  predicted  (optimal)  infiltration  routes 
between  various  start  and  specific  end  points  within  a  particular  area 
of  study. 

The  ability  to  model  real  infiltration  routes  rests  heavily  on 
the  assumption  that  there  is  a  rationale  underlying  the  choice  of  such 
routes,  the  criterion  (or  criteria)  for  which  can  be  discovered  by  the 
modeler.  The  criterion  usually  involves  the  minimization  of  some  costs, 
and  where  these  costs  can  be  expressed  quantitatively,  the  techniques 
of  dynamic  programming  can  be  applied  to  determine  minimum-cost  paths. 

An  inefficient  infiltration  route  may  incur  many  additional  costs 
for  the  infiltrator,  such  as  having  to  spend  more  time  on  the  route, 
carry  more  food  and  water,  endure  more  physical  and  psychological  stress 
and  undergo  greater  exposure  to  a  variety  of  risks.  The  model  presented 
here  expresses  many  of  these  costs  in  terms  of  caloric  requirements, 
and  the  resulting  predicted  routes  represent  a  minimization  of  this 
metric. 

Although  we  experimented  with  cost  functions  that  included  both  a 
linear  and  a  quadratic  dependence  on  grade  in  order  to  allow  for  the 
lower  energy  requirement  of  small  negative  grades,  the  results  were  not 
entirely  satisfactory  because  of  uncertainty  as  to  the  relative  weights 
of  the  linear  and  quadratic  terms.  The  usefulness  of  a  linear  term  in 
the  cost  function  is  further  nullified  because  between  the  same  start 
and  end  points  a  linear  term  contributes  the  same  cost  to  all  routes, 
since  the  total  cost  due  to  this  term  depends  only  on  the  altitude  dif¬ 
ference  between  the  start  and  end  points. 

Since  infiltration  routes  are  traveled  in  both  directions  (because 
of  personnel  rotation  and  transportation  of  supplies  from  either  direc¬ 
tion)  ,  and  since  an  antisymmetric  linear  term  does  not  contribute  to 
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route  differentiation,  the  final  form  of  the  cost  function  employed  in 
the  model  is  symmetric  in  the  cost  associated  with  the  altitude  change 
along  a  predicted  route  (i.e.,  the  cost  associated  with  that  part  of 
the  route  which  goes  downhill  is  the  same  as  the  cost  associated  with 
that  part  going  uphill,  traveling  in  the  same  direction  for  equal  dis¬ 
tances).  This  more  useful  form  is 

Cost  =  1  +  (a  |grade|)^. 

Values  for  p  of  1  and  2  were  investigated.  Values  investigated  for  a 
were  10,  20,  and  40. 

The  model  presently  uses  only  terrain  data  in  its  cost  function 
but  is  easily  expanded  to  include  other  data  as  they  become  available. 
Among  the  necessary  inputs  to  the  model  are  altitude  readings  and  the 
specification  of  the  terminal  point,  all  within  the  area  under  study. 

A  two-dimensional  interpolation  function  may  be  employed  to  interpolate 
for  altitude  values  between  the  input  grid  data  points,  resulting  in  a 
finer  network  representation  of  the  study  area.  The  minimum-cost  paths 
are  computed  from  all  included  points. 

Since  the  optimal  path  from  any  point  is  unique  and  independent  of 
the  path  to  that  point,  it  is  sufficient  for  the  computer  to  print  out 
a  table  that  specifies  to  which  of  the  eight  adjacent  points  one  pro¬ 
ceeds  from  any  point  on  the  grid  en  route  to  the  terminal.  Also  ob¬ 
tained  at  each  point  is  a  cost  measure  representing  the  total  computed 
cost  from  each  point  to  the  terminal. 

The  model  has  been  applied  to  a  rectangular-shaped  area,  210  x  220 
km.  This  grid  covers  part  of  northern  South  Vietnam  and  adjacent  Laos. 
The  road  network  in  Laos  was  used  as  the  location  for  start  points  for 
determining  infiltration  routes  which  terminated  in  various  assumed 
base-camp  areas  in  South  Vietnam.  The  results  obtained  from  this  first 
application  resembled  many  characteristics  of  observed  infiltration 
routes  and  offered  an  explanation  for  some  route  characteristics  not 
previously  understood. 

The  model  is  potentially  widely  applicable.  The  main  areas  of  ap¬ 
plication  are  as  follows: 

1.  Infiltration  route  characteristics ♦  It  has  never  been  seriously 
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suggested  that  infiltration  routes  are  chosen  at  random.  Traditional 
communication  routes  almost  always  follow  easy  rather  than  difficult 
terrain.  The  model  provides  a  systematic  approach  for  investigating 
the  criteria  and  degree  of  randomness  involved  in  route  determination, 
provided  some  infiltration  routes  are  known  to  the  investigator.  The 
model  can  be  useful  in  indicating  probable  infiltration  routes  and  pat¬ 
terns;  identifying  likely  alternative  routes  (for  example,  after  the 
user  has  blocked  optimal  routes);  and  estimating  logistics  and  time  re¬ 
quirements  for  infiltration. 

If  there  were  some  systematic  departure  of  the  real  route  from  the 
"terrain  optimal"  routes  in  some  region,  an  investigation  would  then  be 
made  as  to  what  conditions  in  the  region  caused  the  routes  to  differ. 
Perhaps  the  existence  of  a  village,  a  traditionally  established  path, 
vegetation  cover,  or  unknown  base  camps  could  explain  the  deviation. 

Once  a  reason  for  the  deviation  is  conjectured,  a  new  cost  function  can 
be  constructed  which  includes  a  weighting  factor  for  the  new  condition(s) 

Route  characteristics  other  than  the  optimizing  criteria  can  be  in¬ 
ferred.  For  example,  it  may  be  possible  to  discover  camps  or  rallying 
points  by  examining  the  deviation  of  real  routes  from  those  predicted. 

Or  it  may  be  possible  to  determine  start  or  end  points  for  the  real 
routes  by  matching  fragmentary  route  information  to  routes  predicted 
using  an  array  of  start  and  end  points. 

2.  Barrier  and  sensor  emplacement .  The  model  provides  a  useful 
indicator  of  the  best  places  for  locating  barriers,  sensors,  or  inter¬ 
diction  campaigns  and  allows  for  the  study  of  cost  sensitivity  to  changes 
in  routes  after  blocking  or  detection. 

3.  Offensive  applications .  The  model  should  be  of  value  to  the 
infiltrator  as  well  as  the  defender,  especially  where  the  infiltrating 
force  has  little  experience. 

The  model  can  also  be  applied  in  the  areas  of  selecting  routes  for 
new  roads  or  for  ground  resupply;  finding  ambush  points;  or  determining 
boundaries  of  political  subdivisions  on  the  basis  of  ease  of  communica¬ 
tion — an  area  which  is  too  broad  for  further  consideration  at  this  time. 
The  model  is  conceptually  simple,  but  it  is  diverse  in  the  parameters 
and  implications  it  encompasses  and  therefore  provides  a  rich  field 
for  further  investigation. 
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I.  INTRODUCTION 


BACKGROUND 

Many  of  the  political  boundaries  separating  long-established  coun¬ 
tries  of  the  world  lie  in  rugged  terrain.  The  barrier  that  such  ter¬ 
rain  imposed  on  communications  was  often  responsible  for  the  cultural 
and  political  differentiation.  Today,  when  warfare  occurs  between  coun¬ 
tries  where  such  barriers  exist,  it  may  maintain  many  primitive  aspects, 
particularly  with  regard  to  logistics  and  communications  across  these 
boundaries,  in  spite  of  modern  technological  developments  provided  by 
new,  larger  powers.  Because  modern  implements  of  warfare  such  as  air¬ 
craft  and  tanks  cannot  operate  effectively  over  rugged  terrain  covered 
with  thick  vegetation,  insurgent  groups  fighting  civil  wars  often  use 
such  terrain  for  their  base  camps  and  hidden  supply  routes.  Therefore, 
transportation  and  communication  through  rough  terrain  on  foot  or  by 
small  animal  or  bicycle  pack  remains  an  important  aspect  of  modern  war¬ 
fare  at  the  lower  end  of  the  technological  spectrum. 

An  initial  model  for  simulating,  analyzing,  and  predicting  infiltra¬ 
tion  routes  has  been  completed.  PREMIR  (Prediction  Model  for  Infiltra¬ 
tion  Routes)  is  primarily  a  mathematical  computer  model  which  by  means 
of  established  analytical  techniques  arrives  at  a  set  of  predicted 
(optimal)  infiltration  routes  between  various  start  and  specific  end 
points  within  a  particular  area  of  study.  Although  the  model  presently 
uses  only  terrain  data  in  its  cost  function,  it  can  be  easily  expanded 
to  include  other  data  as  they  become  available. 

Travel  in  rough  terrain  is  more  difficult  than  is  generally  thought 

by  those  who  have  little  experience  in  such  travel.  As  a  rule  of  thumb, 

it  requires  10  times  more  energy  to  ascend  1  foot  in  altitude  than  to 

* 

cover  1  foot  on  horizontal  terrain.  This  difference  in  energy  require¬ 
ments  becomes  even  greater  when  a  heavily  laden  bicycle  is  being  pushed 
up  an  incline  rather  than  over  flat  land.  Consequently,  all  possible 
infiltration  routes  through  rugged  terrain  are  not  equally  suitable 
and,  more  important,  are  not  equally  costly. 


Or  twice  as  much  energy  to  ascend  a  10  percent  grade  as  to  cover 
a  horizontal  path. 
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Of  course,  terrain  alone  is  not  the  only  important  factor  in  choos¬ 
ing  infiltration  routes.  The  availability  of  vegetation  to  hide  the 
routes,  the  availability  of  food  and  water,  the  avoidance  of  unfriendly 
population,  and  the  availability  of  friendly  population  will  also  alter 
an  infiltration  route  based  upon  terrain  alone. 

The  idea  of  an  infiltration  route  assumes  that  there  is  a  definite 
start  point  and  end  point.  However,  in  the  course  of  a  conflict,  these 
points  may  change  and  the  corresponding  routes  may  change,  so  that  one 
must  be  careful  in  assuming  that  such  routes  have  a  permanent  character 
or  that  the  start  points  and  end  points  will  be  known  in  advance.  Never¬ 
theless,  in  order  to  develop  a  model  for  predicting  infiltration  routes, 
it  is  necessary  to  make  certain  idealizing  assumptions,  and  one  of  the 
first  assumptions  usually  employed  is  that  terminal  points  or  initial 
points  are  known.  Later  on  we  will  be  able  to  relax  even  this  assump¬ 
tion. 

The  ability  to  model  real  infiltration  routes  rests  on  the  assump¬ 
tion  that  there  is  a  rationale  behind  choosing  such  routes,  the  criteria 
for  which  can  be  discovered  by  the  modeler.  In  other  words,  infiltra¬ 
tion  routes  are  not  chosen  completely  at  random  but  are  chosen  to  mini¬ 
mize  some  costs,  whether  they  be  physical,  economic,  or  psychological. 

In  the  model  we  have  developed,  it  is  assumed  that  these  cost  factors 
can  be  measured  by  a  common  metric.  For  example,  one  might  be  willing 
to  walk  10  miles  out  of  his  way  to  avoid  scaling  a  1-mile-high  peak, 
or  to  enable  50  miles  of  his  path  to  lie  beneath  a  jungle  canopy  rather 
than  in  the  open,  or  to  avoid  passing  within  1  mile  of  an  unfriendly 
village. 

With  such  cost  factors  based  on  a  common  metric  it  is  possible  to 
utilize  a  dynamic-programming  algorithm  which  will  determine  the  minimum- 
cost  path  between  two  points.  These  predicted  paths  can  be  used  in  many 
ways.  For  example,  one  can  assume  a  start  point  such  as  a  road  in  Laos 
and  an  end  point  such  as  a  base  camp  in  South  Vietnam  and  find  the 
minimum-cost  infiltration  route  based  on  the  assumed  form  of  the  cost 
function.  If  segments  of  an  infiltration  route  are  already  known  through 
aerial  reconnaissance  or  prisoner  interrogations,  these  segments  can  be 
compared  to  the  prediction,  and  differences  can  be  resolved  by  trying 
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new  cost  functions  or  by  trying  new  start  and  end  points.  In  this  man¬ 
ner,  we  have  a  systematic  way  of  investigating  the  rationale  used  by 
the  infiltrator  in  choosing  his  infiltration  routes. 

The  model  can  also  be  used  to  locate  optimal  points  for  construct¬ 
ing  infiltration  barriers  or  for  interdicting  infiltration.  Barriers 
can  be  introduced  into  the  model  mathematically  by  assigning  high  costs 
to  certain  points  across  an  infiltration  route  and  observing  the  new 
infiltration  routes  predicted  by  the  model  which  circumvent  these  bar¬ 
riers.  A  comparison  of  the  costs  to  the  infiltrator  before  and  after 
barrier  installation  gives  a  measure  of  the  effectiveness  of  the  spe¬ 
cific  barrier  installation.  One  would  expect  barriers  to  be  most  effec¬ 
tive  at  choke  points,  that  is,  regions  where  many  infiltration  routes 
with  different  start  and  end  points  pass  close  to  one  another  because 
of  the  high  cost  associated  with  alternate  areas.  These  choke  points 
show  up  quite  clearly  in  some  of  the  model  calculations  that  have  been 
performed. 

The  area  from  which  terrain  data  were  gathered  included  the  north¬ 
ern  portion  of  South  Vietnam  and  the  adjacent  areas  of  Laos.  The  road 
network  in  Laos  was  used  as  the  location  for  start  points  for  determin¬ 
ing  infiltration  routes  which  terminated  in  various  base-camp  areas  as¬ 
sumed  in  Vietnam.  The  results  obtained  resembled  many  characteristics 
of  observed  infiltration  routes  and  offered  an  explanation  for  some  of 
the  route  characteristics  not  previously  understood. 

FACTORS  INFLUENCING  CHOICE  OF  INFILTRATION  ROUTES 

An  inefficient  infiltration  route  may  incur  many  costs  associated 
with  factors  such  as  the  following:  The  infiltrator  may  have  to  spend 
more  time  en  route;  carry  more  food  or  water;  endure  more  physical  and 
psychological  stress;  risk  greater  exposure  to  disease  as  a  result  of 
insect  or  animal  bites;  risk  greater  chance  of  discovery  by  the  enemy; 
and  risk  greater  chance  of  losing  his  way  and  not  reaching  his  destina¬ 
tion. 

The  food  and  water  requirement  may  be  the  most  convenient  means 
of  measuring  the  overall  cost  factors  associated  with  physical  stress. 

In  many  parts  of  the  world,  springs  and  streams  are  readily  available, 
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so  that  water  need  not  be  carried  by  the  infiltrator  except  in  small 
amounts.  Food,  on  the  other  hand,  is  not  often  available  in  rugged 
terrain  and  is  usually  carried  by  the  infiltrator;  in  Southeast  Asia 
rice  is  most  commonly  carried.  In  Indochina,  the  infiltrator,  who  may 
be  on  foot  for  as  long  as  3  to  6  months,  generally  cannot  carry  with 
him  all  the  food  that  will  be  required  for  such  a  trip.  Therefore, 
logistics  bases  and  rest  areas  will  usually  be  set  up  along  the  infil¬ 
tration  route,  giving  the  route  a  certain  characteristic  of  permanence, 
since  it  is  not  easy  to  change  the  logistics  system  frequently. 

The  cost  of  this  logistics  system  itself  will  be  proportional  to 
the  rate  at  which  food  is  consumed  by  individual  infiltrators  and  to  the 
total  number  of  infiltrators.  Since  the  food  will  usually  be  brought 
into  the  supply  areas  by  some  primitive  methods,  the  rough  terrain  im¬ 
poses  an  additional  burden  on  the  logistics  system  as  well  as  on  the 
infiltrator.  In  fact,  the  choice  of  infiltration  routes  may  be  con¬ 
strained  more  by  the  need  for  transportation  of  supplies  to  the  rest 
areas  along  the  trail  than  by  the  stress  placed  on  the  infiltrator  as 
an  individual. 

Since  supplies  may  have  to  move  in  a  direction  opposite  to  that 
of  infiltration,  such  as  from  the  rice  paddies  of  Vietnam  up  into  the 
highlands,  we  cannot  assume  that  costs  should  be  assigned  only  to  that 
part  of  the  infiltration  route  which  is  downhill  for  the  infiltrator. 

For  this  reason,  in  our  final  form  of  the  cost  function,  we  have  as¬ 
sumed  symmetry  in  the  cost  associated  with  altitude  change  along  an 
infiltration  route;  that  is,  the  cost  of  that  part  of  the  route  which 
goes  downhill  is  as  great  as  the  cost  of  that  part  going  uphill.  This 
form  of  the  function,  however,  is  not  final  and  may  be  changed  if  evi¬ 
dence  indicates  it  should  be  otherwise. 

For  the  reasons  discussed,  we  have  assumed  that  the  amount  of  food 
required  per  infiltrator  will  be  an  important  factor  in  determining  the 
infiltration  routes  employed.  The  rate  at  which  food  is  consumed  will 
depend  upon  a  number  of  variables  which  we  will  discuss. 

We  assume  that  food  is  consumed  primarily  for  calories,  other  nu¬ 
trients  being  obtained  in  the  process  with  little  additional-cost.  The 
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calories  required  per  mile  of  travel  will  depend  upon  the  speed  of 
travel.  An  infiltrator  running  along  flat  land  consumes  more  calories 
per  mile  than  one  traveling  more  slowly  at  a  comfortable  speed  of  2 
or  3  miles  per  hour.  Traveling  too  slowly  also  requires  greater  amounts 
of  food  per  mile.  Traveling  up  a  slope  requires  more  calories  than 
traversing  flat  land.  Data  provided  later  in  this  Report  indicate  that 
10  times  as  many  calories  are  required  for  every  foot  of  altitude  gained 
as  are  required  for  1  foot  of  horizontal  travel. 

The  roughness  of  the  terrain  will  also  be  a  factor  in  determining 
the  difficulty  of  the  path  and  hence  the  amount  of  food  required.  It 
is  easier  to  walk  a  smooth  path  on  flat  land  than  it  is  to  walk  along 
a  ridgeline  or  along  the  bottom  of  a  narrow  valley,  even  though  no  net 
change  in  altitude  is  involved.  Furthermore,  ridgelines  and  valleys 
are  often  easier  to  follow  than  paths  following  the  contour  lines  on 
steep  hillsides,  since  one  must  exert  additional  energy  to  obtain  solid 
footing  and  to  avoid  slipping  down  the  hill.  Even  if  the  paths  are 
well  prepared,  erosion  will  often  remove  level  surfaces  from  sections 
of  the  path. 

Temperature,  humidity,  and  altitude  have  also  been  found  to  affect 
calorie  requirements.  When  temperatures  are  very  low,  additional  food 
is  required  to  maintain  body  temperature;  when  temperatures  are  very 
high,  metabolic  efficiency  is  reduced  and  additional  calories  are  re¬ 
quired  to  do  the  same  work  done  in  moderate  temperatures.  High  humid¬ 
ity  causes  reduction  in  efficiency  of  heat  loss  from  the  body  and  a 
consequent  rise  in  body  temperature  to  an  inefficient  level.  During 
conditions  of  very  low  humidity,  body-water  loss  is  increased,  and  addi¬ 
tional  amounts  of  water  may  have  to  be  carried  to  replace  this  loss. 

The  metabolic  efficiency  of  a  person  not  adapted  to  high  altitudes  may 
be  reduced  when  he  must  function  in  such  areas. 

All  of  these  factors  are  augmented  if  additional  weight  must  be 
carried  by  the  infiltrator.  For  example,  if  there  is  a  long  distance 
between  logistics  bases,  a  considerable  amount  of  additional  food  and 
water  may  have  to  be  carried,  and  in  turn  more  food  and  water  is  re¬ 
quired  to  carry  the  additional  weight,  so  that  for  long  distances  the 
weight  to  be  carried  increases  more  than  linearly  with  the  distance 
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traversed.  The  problem  becomes  even  more  severe  if  the  weight  carried 
is  increased  by  weapons  and  ammunition. 

It  would  appear  that  establishing  infiltration  routes  through  the 
jungle  is  less  of  a  problem  than  was  first  thought.  Most  jungles  are  , 
heavily  populated  with  animals,  which  have  developed  trails  for  their 
own  movement  through  the  jungle.  Some  of  these  trails  cover  consider¬ 
able  distance — often  in  favorable  directions  for  the  infiltrators — as 
seasonal  factors  require  changes  in  feeding  habits.  Furthermore,  in 
situations  such  as  that  in  South  Vietnam,  where  infiltration  routes 
have  been  developed  over  years  and  tens  of  thousands  of  personnel  have 
traversed  these  routes,  the  initial  cost  of  clearing  routes  that  are 
optimal  on  the  basis  of  other  considerations  is  small. 

The  jungle  can  actually  be  a  preferred  area  for  infiltration  routes 
because  of  the  cover  it  provides  against  detection  from  the  air.  How¬ 
ever,  in  Vietnam  where  jungle  covers  are  plentiful,  many  infiltration 
routes  are  found  to  pass  through  open  areas.  This  may  imply  that  the 
need  for  air  cover  is  not  an  important  determinant  in  the  selection  of 
infiltration  routes. 

The  presence  of  friendly  or  unfriendly  villages  may  alter  the 
course  chosen  for  an  infiltration  route.  Friendly  villages  can  serve 
as  rest  and  resupply  areas  for  the  infiltrator  and  as  a  source  of  food 
and  water.  Establishment  of  infiltration  routes  near  unfriendly  vil¬ 
lages  can  result  in  discovery  of  these  routes  and  would  thus  be  avoided 
by  those  establishing  them.  However,  it  is  probably  important  to  con¬ 
sider  that  in  the  course  of  civil  war,  the  loyalties  of  various  vil¬ 
lages  may  shift  from  one  side  to  the  other,  and  even  a  friendly  village 
may  have  an  informer;  thus  it  would  perhaps  be  in  the  best  interest 
of  the  infiltrator  to  establish  his  important  routes  far  from  any  popu¬ 
lation  centers.  Such  a  choice  may  be  very  difficult  in  some  cases, 
however,  since  villages  in  mountainous  areas  will  probably  be  connected 
by  long-established  communication  routes  which  may  very  well  be  the  best 
and  most  easily  passable  in  the  area,  particularly  if  the  terrain  is 
rugged. 

A  river  can  be  either  a  barrier  or  an  aid  to  infiltration,  depend¬ 
ing  upon  its  type  of  flow  and  direction  relative  to  the  direction  of 
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inf iltration.  In  the  lowlands,  rivers  tend  to  be  broad  and  slow-moving, 
and  they  leave  flat  valleys  which  are  easily  followed.  However,  in 
much  of  the  world  and  particularly  in  Southeast  Asia,  these  river  val¬ 
leys  (because  of  their  fertility  and  the  ease  of  communication  offered 
by  the  river)  are  often  heavily  populated,  so  that  use  of  the  river  by 
the  infiltrator  involves  passing  through  populated  areas,  with  the  im¬ 
plications  described  above. 

In  mountainous  country,  rivers  tend  to  be  narrow  and  swift,  and  the 
passages  cut  through  the  mountains  are  often  narrow,  with  no  adjacent 
footpaths  along  the  shore.  In  Southeast  Asia,  the  rivers  do  not  seem 
to  form  difficult  barriers  for  most  infiltration  routes,  since  they  are 
easily  crossed  by  rope  suspension  bridges  or  by  improvised  ferry  ser¬ 
vice.  Furthermore,  a  swift  river  may  be  an  advantage  in  that  it  often 
provides  a  good  source  of  fresh  water. 

Seasonal  factors  may  play  an  important  role  in  route  determination. 
Routes  which  are  very  good  during  one  season  of  the  year  may  become  im¬ 
passable  during  others.  During  the  rainy  season  in  Southeast  Asia,  for 
example,  some  infiltration  routes  are  completely  covered  by  flood  waters; 
others  may  become  so  thick  with  mud  that  it  is  impossible  to  move  bi¬ 
cycles  along  them,  and  walking  requires  additional  energy  consumption. 
Rain  may  also  seriously  interfere  with  the  transportation  of  food  such 
as  rice,  which  deteriorates  in  moisture.  Serious  hazards  are  presented 
in  desert  areas  during  the  infrequent  but  heavy  rains.  Because  of  the 
bareness  of  the  hillsides,  water  collects  very  rapidly  in  arroyos  and 
may  move  as  a  wall  several  feet  high  down  narrow  gulleys  at  a  high  rate 
of  speed.  Passage  through  narrow  canyons  in  the  desert  during  the  rainy 
season  can  therefore  be  very  hazardous,  and  infiltration  routes  at  higher 
altitudes  might  be  preferred. 

During  the  summer  season  passage  through  desert  areas  may  become 
difficult  due  to  the  extreme  heat,  requiring  the  infiltrators  to  move 
only  at  night.  Traffic  at  night,  however,  must  be  accomplished  in  com¬ 
plete  darkness,  since  infiltrators  are  readily  spotted  at  great  dis¬ 
tances  if  they  are  using  any  light  to  illuminate  their  path.  Conse¬ 
quently,  infiltration  routes  for  night  travel  through  desert  areas  must 
be  chosen  with  more  caution.  A  greater  cost  would  therefore  be  assigned 
to  difficult  mountainous  routes. 
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The  unavailability  of  water  is  always  a  serious  problem  in  the 
desert,  especially  during  the  summer  season,  since  considerable  quan¬ 
tities  of  water  may  be  required  to  replace  the  high  water  loss  of  the 
body.  Additional  water  must  therefore  be  carried  and  thus  a  greater 
cost  is  incurred.  One  would  expect  infiltration  routes  to  be  altered 
to  take  advantage  of  any  water  availability  along  the  way. 

In  winter  different  kinds  of  problems  are  presented  to  the  infil¬ 
trator.  Extreme  cold  may  require  that  he  build  fires  to  prepare  food 
and  to  keep  warm  during  the  night.  In  such  cases,  covered  terrain  will 
be  more  desirable  in  order  that  the  fires  will  not  be  readily  detected. 
The  cover  may  be  thick  vegetation  or  caves  and  rocky  shelters.  Heavy 
snowfall  may  make  mountain  canyons  impassable  for  certain  kinds  of 
traffic.  On  the  other  hand,  in  ski  country,  some  passages  may  be  ren¬ 
dered  considerably  easier  by  the  presence  of  snow,  troops  being  able 
to  cover  many  more  miles  per  day  on  skis  than  they  could  on  foot  in  the 
absence  of  snow. 

Of  the  seasonal  factors  discussed  above,  only  two  would  be  signif¬ 
icant  in  affecting  choice  of  infiltration  routes  in  the  Indochina  area 
chosen  for  application  of  this  model:  rains  which  flood  infiltration 
routes  and  deteriorate  foodstuffs,  and  winter  cold  which  requires  the 
infiltrator  to  build  fires  for  warmth  and  food  preparation. 
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II,  DATA  COLLECTION  AND  INTERPOLATION 

In  order  to  apply  the  model  to  a  real  infiltration  problem,  we 
chose  a  rectangular  area  in  Indochina,  210  *  220  km.  This  area  in¬ 
cludes  part  of  northern  South  Vietnam  and  adjacent  Laos  lying  between 
15°  and  17°  latitude  and  106°  and  108°  longitude  (Fig.  1). 

The  dynamic-programming  technique  to  be  used  for  minimum-cost  path 
determination  requires  that  a  finite  set  of  nodal  points  be  designated 
and  the  cost  from  traveling  from  any  point  to  an  adjacent  point  be  spec¬ 
ified.  The  number  of  nodal  points  is  limited  by  computer  capacity  and 
costs  of  data  collection  and  computation. 

For  convenience,  we  chose  a  rectangular  grid  of  nodal  points.  The 
distance  between  nearest  points  was  set  at  2.5  km  because  of  computer 
capacity.  Because  of  the  time  required  to  collect  altitude  data  over 
a  grid  of  approximately  8000  points,  it  was  decided  to  collect  data  at 
a  smaller  number  of  points  and  interpolate. 

Several  techniques  were  tried  for  collecting  the  altitude  data  in 
digital  form.  A  graphical  input  device,  the  Rand  Tablet,  was  used  to 
trace  contour  lines  from  maps  of  scale  1:250,000,  altitudes  at  grid 
points  being  determined  by  interpolation  from  altitudes  where  contour 
lines  crossed  grid  lines.  It  was  found,  however,  that  this  tracing 
technique  was  very  time-consuming,  and  the  required  programming  for  in¬ 
terpolation  was  more  complex  than  first  appeared. 

The  simplest  way  to  collect  altitude  data  was  the  one  finally  used. 

A  set  of  plastic  relief  maps  of  the  area  chosen  at  a  scale  of  1:250,000 
was  employed.  These  maps  had  altitude  contours  printed  at  100-m  inter¬ 
vals.  Each  500-m  contour  line  was  labeled.  The  raised  relief  feature 
of  the  map  allowed  us  to  determine  quickly  whether  a  given  contour  line 
was  higher  or  lower  in  altitude  than  a  labeled  line.  Also  printed  on 
the  maps  were  the  standard  10-km  grid  lines.  The  points  at  the  inter¬ 
sections  of  the  grid  lines  were  chosen  for  altitude  readings.  Altitude 
at  grid  points  was  estimated  to  the  nearest  10  m,  by  interpolating  by  eye, 
using  the  distance  from  the  grid  point  to  adjacent  100-m  contour  lines. 

Because  there  were  506  grid  points  at  which  altitude  data  were  ob¬ 
tained  and  more  than  8000  points  at  which  altitude  data  were  required, 
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Fig.l  —  Northern  region  of  South  Vietnam  and  adjacent  Laos 
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several  interpolation  schemes  were  investigated.  This  investigation 

; k 

was  conducted  by  using  a  JOSS  program  which  allowed  quick  implementa¬ 
tion  of  an  interpolation  technique  and  provided  a  map  printout  on  which 
contour  lines  could  be  drawn  and  compared  to  those  on  a  geodetic  map. 

The  final  interpolation  technique  chosen  because  of  its  resemblance 
to  the  real  map  is  analogous  to  the  SYMAP^"^  technique  developed  at  Har¬ 
vard  for  interpolation  of  variables  over  two-dimensional  surfaces.  Each 
point  to  be  interpolated  is  given  an  altitude  which  is  a  weighted  aver¬ 
age  of  the  altitudes  of  the  12  adjacent  altitudes  on  the  primary  10-km 
grid.  The  weighting  function  is  the  inverse  square  of  the  distance. 

In  other  words,  the  altitude  of  the  point,  P,  under  consideration  is 
proportional  to  the  sum  of  the  altitudes  at  each  adjacent  grid  point 
divided  by  the  square  of  the  distance  between  P  and  the  grid  point. 

The  weighted  sum  is  then  divided  by  the  sum  of  the  reciprocals  of  the 
squares  of  the  distances  in  order  that  the  weighting  factors  add  up  to  1. 

This  interpolation  technique  has  a  certain  intuitive  appeal.  If 
a  cone-shaped  mountain  were  to  erode  to  a  cone  with  a  larger  base,  the 
altitude  at  a  point  halfway  from  the  center  of  the  cone  would  decrease 
as  the  inverse  square  of  the  distance  of  this  point  from  the  center. 

One  disadvantage  of  this  interpolation  technique  is  that  no  alti¬ 
tude  at  an  interpolated  point  can  be  a  local  maximum  (or  minimum)  point. 
That  is,  no  mountain  peak  interior  to  the  grid  can  be  predicted  by  this 
technique.  Interpolation  techniques  that  extrapolated  altitude  trends 
to  correct  this  deficiency  were  found  to  yield  unrealistic  results  be¬ 
cause  the  terrain  has  considerable  random  structure  at  scales  smaller 
than  the  10-km  grid  size. 

The  failure  of  the  technique  to  predict  local  minima  or  basins  was 
not  viewed  as  a  handicap  because  of  the  relative  absence  of  these  fea¬ 
tures  in  the  terrain  being  examined.  The  ability  of  the  technique  to 
structure  ridges  and  valleys  accurately  is  probably  very  poor,  however, 

_ 

The  JOSS  program  is  given  in  Appendix  A.  JOSS  is  the  trademark 
and  service  mark  of  The  Rand  Corporation  for  its  computer  program  and 
services  using  that  program.  The  JOSS  system  is  an  on-line,  time-shared 
system  developed  at  Rand. 
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and  this  is  a  major  concern  for  infiltration-route  prediction  because 
of  their  importance.  This  can  be  compensated  for  to  a  certain  extent 
by  manipulation  of  the  cost  function,  which  we  will  discuss  later. 

As  an  example  of  the  ability  of  the  inverse-square  interpolation 
technique  to  reproduce  terrain,  we  will  consider  a  40  x  40  km  area  of 
South  Vietnam  containing  the  A  Shau  valley.  Part  (a)  of  Fig.  2  is  a 
reproduction  of  the  300-m  contour  lines  drawn  from  a  map  of  scale 
1:1,000,000  and  enlarged.  The  elevation  in  meters  is  shown  by  the  ac¬ 
companying  legend.  The  ,!T"-shaped  region  above  center  between  300  and 
600  m  elevation  is  a  segment  of  the  A  Shau  valley.  The  low  area  (less 
than  300  m  elevation)  in  the  upper  right-hand  corner  of  the  figure  leads 
to  the  coastal  plain. 

Part  (b)  of  Fig.  2  shows  a  JOSS  printout  of  the  altitudes  covering 
* 

the  same  area. 

The  single  digit  in  every  other  column  represents  the  altitude 
rounded  to  100  m.  For  example,  the  5  in  the  lower  left-hand  corner  in¬ 
dicates  the  altitude  at  this  point  is  between  500  and  600  m.  We  drop 
the  first  digit  for  altitudes  above  1000  m;  for  example,  4,  the  19th 
digit  in  the  bottom  row,  represents  1400  to  1500  m. 

Altitudes  were  read  into  the  program  at  the  10-km  UTM  (Universal 
Transverse  Mercator)  grid  intersections  in  Fig.  2(a)  as  well  as  at  the 
24  intersections  surrounding  the  figure.  Altitudes  were  interpolated 
on  a  1-km  grid  at  all  other  points,  using  the  inverse-square  technique. 

By  comparing  (a)  and  (b)  of  Fig.  2,  we  can  see  that  the  interpola¬ 
tion  program  follows  the  gross  terrain  characteristics  rather  well.  The 
high-altitude  regions  at  the  bottom  of  part  (a)  are  reproduced  as  a  re¬ 
sult  of  their  falling  on  grid  intersections.  The  low-altitude  region 
in  the  upper  right-hand  corner  is  well  reproduced.  The  valley  areas, 
however,  are  not  reproduced  as  well  as  was  expected. 

The  altitudes  on  the  10-km  grid  were  obtained  from  maps  of  larger 
scale  (1:250,000)  than  that  of  Fig.  2(a).  It  has  been  observed  that 
these  maps  do  not  always  agree  on  the  altitude  at  a  grid  point,  and 

Only  the  digits  in  the  odd-numbered  columns  are  significant.  The 
I’s  in  the  even-numbered  columns  are  inserted  only  for  formatting  pur¬ 
poses. 
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(a )  Map  contours 


(b)  Interpolation  results 
Fig. 2 — Interpolation 
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this  will  account  for  some  of  the  discrepancy  between  the  interpolated 
map  and  the  maps  used  for  illustration  in  this  Report,  It  has  also  been 
found  that  maps  of  the  same  scale  do  not  always  agree,  and  this  is  due 
in  part  to  uncertainties  in  altitude  of  some  of  the  regions  of  Vietnam 
when  the  maps  were  prepared.  Because  terrain  data  used  in  this  study 
are  of  limited  accuracy,  details  of  the  routes  predicted  by  the  model 
are  not  to  be  taken  very  seriously.  The  data  used  are  considered  suf¬ 
ficient  to  illustrate  the  value  and  applications  of  this  approach  but 
are  insufficient  for  strategic  or  tactical  use. 

Although  our  discussion  has  been  concerned  primarily  with  terrain 
data,  other  kinds  of  data  can  be  quantified  and  interpolated  in  the  same 
manner.  For  example,  at  each  grid  point  an  integer  can  be  assigned  in¬ 
dicating  the  amount  of  vegetative  cover  for  the  infiltration  route  in 
the  area.  Heavy  jungle  can  be  assigned  a  value,  1;  light  cover  with 
occasional  trees  or  brush  can  be  assigned  a  value,  2;  no  cover  can  be 
assigned  a  value,  3.  And  these  values  can  be  interpolated  on  the  grid 
in  the  same  manner  as  the  altitude  data.  The  interpolated  values  can 
then  be  used  in  the  route  cost  function  as  described  later. 

Population  data  can  be  handled  in  a  slightly  different  fashion. 

The  map  coordinates  and  population  of  each  village  can  be  read  into  the 
computer  program.  A  number  can  be  assigned  to  each  interpolation  point 
that  gives  a  measure  of  population  density  by  adding  the  populations  of 
each  village  within  a  given  radius  of  the  point  (10  km,  for  example) 
weighted  by  the  inverse  square  of  the  distance  from  that  village  to  the 
point.  The  resulting  number  can  be  used  as  a  population  cost  factor 
when  routes  are  chosen  to  avoid  (or  seek)  population.  A  positive  or 
negative  factor  could  be  used  for  each  village,  depending  on  whether  it 
is  regarded  as  unfriendly  or  friendly.  Similarly,  one  can  treat  other 
geographically  distributed  factors  that  may  be  important  to  choice  of 
infiltration  route. 


Ill .  THE  COST  EQUATION' 

The  central  hypothesis  of  our  model  is  that  infiltration  routes 
are  chosen  to  minimize  some  cost,  and  that  this  cost  accumulates  incre¬ 
mentally  along  the  route.  In  other  words,  the  routes  are  solutions  to 
minimum-path  problems. 

In  order  to  solve  a  minimum-path  problem  by  a  dynamic-programming 
algorithm,  it  is  necessary  to  specify  a  cost  equation  or  weighting 
function,  giving  the  incremental  costs  along  arcs  in  the  grid.  We 
might  expect  that  in  the  real  world  such  a  cost  equation  would  be  a 
complex  function  of  terrain,  vegetation  characteristics,  population 
distribution,  and  other  factors  already  discussed.  Nevertheless,  be¬ 
cause  we  felt  terrain  was  the  primary  factor  determining  infiltration 
routes  in  mountainous  areas,  we  chose  to  examine  the  route  implications 
of  a  cost  equation  based  on  terrain  only.  This  approach  does  not  ex¬ 
clude  other  cost  factors  from  consideration,  since  they  can  be  added 
to  the  cost  equation  once  the  influence  of  terrain  is  understood. 

In  the  cost  equation,  we  chose  to  include  three  terrain  factors: 

o  The  distance  between  adjacent  nodes 
o  The  altitude  difference  between  adjacent  nodes 
o  The  roughness  of  the  terrain  at  each  node 

We  had  a  choice  of  different  ways  to  account  for  the  cost  of  these  fac¬ 
tors  : 

o  The  physical  stress  on  the  infiltrator 

o  The  food  and  water  requirements  of  the  infiltrator  imposed  by 
the  terrain 

o  The  time  required  to  traverse  the  terrain  on  foot 
o  The  costs  to  the  logistics  system  imposed  by  the  terrain 

In  this  case,  we  chose  to  use  as  the  cost  metric  the  number  of  calories 

required  by  the  infiltrator  to  traverse  the  terrain  at  an  efficient 
& 

rate  of  speed.  It  was  felt  that  the  total  stress  on  the  individual 


A  rate  of  speed  at  which  the  number  of  calories  required  per  mile 
of  travel  is  at  a  minimum. 
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would  be  proportional  to  the  work  done  in  traveling.  All  other  things 

being  equal,  the  food  and  water  and  hence  logistics  requirements  would 

be  proportional  to  the  work  done.  Finally,  since  there  is  a  maximum 

rate  at  which  each  individual  can  expend  energy  on  a  continuous  basis 

without  succumbing  to  exhaustion,  the  time  required  for  a  given  route 

should  be  roughly  proportional  to  the  total  energy  expended. 

& 

Table  1  shows  the  rate  of  expenditure  of  energy  for  various  trav¬ 
eling  conditions.  As  would  be  expected,  these  rates  are  dependent  on 

* 

the  speed  of  travel.  This  can  be  seen  more  clearly  in  (a)  of  Fig.  3, 
where  a  line  of  45-deg  slope  represents  a  rate  of  energy  expenditure 
per  unit  distance.  The  most  efficient  walking  speed  is  approximately 
3.5  mi/hr;  at  this  speed,  energy  is  expended  at  a  rate  of  1/20  kcal  per 

•kit 

meter  of  distance  traveled. 

These  data  were  obtained  for  men  of  various  physical  conditions 

traveling  on  flat,  hard  surfaces.  The  effect  of  traveling  uphill  can 

* 

be  inferred  from  treadmill  data,  shown  in  (b)  of  Fig.  3.  It  can  be 
seen  in  this  figure  that  traveling  up  a  10  percent  grade  at  3.5  mi/hr 
requires  nearly  twice  the  energy  expenditure  per  mile  (or  per  hour)  as 
traveling  on  level  ground.  It  is  also  suggested  in  the  running  curves 
that  the  energy  expenditure  (per  mile)  increases  faster  than  linearly 
with  grade.  We  have  examined  both  linear  and  quadratic  dependence  on 
grade,  as  discussed  below. 

A  treadmill  is  a  kind  of  idealized  condition  for  measuring  energy 
expenditure,  since  little  energy  is  required  to  maintain  balance  or 
to  compensate  for  uneven  terrain.  As  the  data  in  Table  1  suggest,  a 
greater  energy  expenditure  is  required  in  stubble  and  plowed  fields 
and  in  rough  country  than  that  shown  in  Fig.  3. 

Most  studies  neglect  to  report  the  energy  requirements  to  descend 
on  a  negative  grade.  There  is  some  evidence,  and  experience  suggests, 
that  a  slight  downhill  grade  requires  less  effort  per  meter  than 

/a\ 

Data  taken  from  Bio  astronautics  Data  Book . 

Since  a  pound  of  rice  contains  approximately  1600  kcal,  we  ob¬ 
tain  a  ’’mileage"  of  32  km  per  pound  of  rice.  This  rate  will  be  smaller 
for  lighter-weight  individuals  carrying  no  pack. 
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Table  1 

ENERGY  EXPENDITURES  OF  ACTIVITIES 


Typical  values  for 

Activity  lb  02/hr  kcal/min  Btu/hr 


Moving  over  rough  terrain  on  foot 
Flat  firm  road 
Grass  path 
Stubble  field 
Deeply  plowed  field 
Steep  45°  slope 
Plowed  field 

Soft  snow,  with  44-lb  load 

Load  carrying 

Walking  on  level 
with  58-lb  load, 
trained  men 

Walking  on  level 
with  67-lb  load, 
trained  men 

Walking  on  level 
with  75-lb  load, 
trained  men 

Walking  up  36%  grade 
with  43-lb  load, 
sedentary  men 

Movement  in  snow 

Skiing  in  loose  snow 
Sled  pulling — low  drag,  hard  snow 
Snowshoeing — bearpaw  type 
Skiing  on  level 

Sled  pulling — low  drag,  medium  snow 

Snowshoeing — trail  type 

Walking,  12-18"  snow,  breakable  crust 

Skiing  on  loose  snow 

Snowshoeing — trail  type 

Skiing  on  loose  snow 


Measured  work  at  different  altitudes 


Bicycle  ergometer  1 

'430  kg-m/min 

430 

(, 

430 

Mountain 

880- 

1037 

kg-m/min 

climbing  ' 

566- 

786 

’ 

393- 

580 

2.5  mph 

0.11-0.19 

2. 8-4. 9 

660-1140 

2.5 

0.12-0.20 

3. 2-5.1 

760-1240 

2.5 

0.16-0.23 

4. 0-6.1 

960-1440 

2.0 

0.19-0.27 

4. 9-6. 9 

1160-1640 

1.5 

0.19-0.27 

4. 9-6. 9 

1160-1640 

3.3 

0.30 

7.8 

1850 

2.5 

0.79 

21.0 

4850 

2 . 1  mph 

0.07 

1.9 

450 

2.7 

0.11 

2.9 

690 

3.4 

0.18 

4.6 

1100 

4.1 

0.32 

8.3 

1980 

2.1 

0.09 

2.3 

550 

2.7 

0.11 

2.9 

690 

3.4 

0.20 

5.1 

1210 

4.1 

0.33 

8.4 

2000 

2.1 

0.10 

2.5 

600 

2.7 

0.13 

3.4 

810 

3.4 

0.20 

5.2 

1240 

4.1 

0.34 

8.6 

2100 

0.5 

0.26 

6.7 

1590 

1.0 

0.47 

12.3 

2910 

1.5 

0.62 

16.0 

3800 

2 . 6  mph 

0.32 

8.1 

1930 

2.2 

0.34 

8.6 

2020 

2.5 

0.34 

8.7 

2070 

3.0 

0.35 

9.0 

2140 

2.0 

0.38 

9.7 

2310 

2.5 

0.40 

10.3 

2460 

2.5 

0.50 

12.7 

3010 

5.2 

0.52 

14.6 

3800 

3.5 

0.59 

14.8 

4200 

8.1 

0.80 

20.6 

4900 

720  mm  Hg 

0.20 

5.1 

1230 

620 

0.19 

4.9 

1170 

520 

0.21 

5.4 

1290 

610  mm  Hg 

0.36-0.43 

9.2-11.0 

2200-2640 

425 

0.30-0.37 

7.7-  9.5 

1840-2260 

370 

0.25-0.41 

6.4-10.5 

1530-2520 

40,000 

30,000 

20,000 


^  10,000 


•§"  5,000 

|  4,000 

<d 

X  3,000 

2,000 


1,000 

700 


~u 

D 

O 


Angle  of  incline  (°) 

(b) 


Speed  (m/sec) 
(a) 


I 


00 

I 


Fig. 3 — Energy  costs  of  foot  travel 
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horizontal  travel.  The  situation  changes,  however,  when  the  downhill 
grade  increases  to  the  point  where  slippage  is  possible.  In  this  case 
considerable  energy  can  be  required  to  maintain  footing. 

Although  we  experimented  with  cost  functions  that  included  both  a 
linear  and  a  quadratic  dependence  on  grade  in  order  to  reflect  the  lower 
energy  requirement  for  small  negative  grades,  the  results  were  not  en¬ 
tirely  satisfactory  because  of  uncertainty  as  to  the  relative  weights 

* 

of  the  linear  and  quadratic  terms.  We  resolved  this  ambiguity  by  argu¬ 
ing  that  infiltration  routes  will  be  traveled  in  both  directions  as  per¬ 
sonnel  are  rotated  and  because  food  is  often  transported  to  camps  in  the 
highlands  from  either  direction.  With  this  argument  we  primarily  exam¬ 
ined  cost  functions  which  depend  only  on  the  absolute  value  of  the  al¬ 
titude  change  between  grid  points  and  which  neglect  the  difference  be¬ 
tween  uphill  and  downhill. 

Some  mountaineering  experts  have  measured  the  energy  required  for 

(3) 

traveling  various  kinds  of  terrain  in  many  kinds  of  weather.  Their 

results  were  obtained  by  use  of  an  oxygen- consumption  meter  strapped 
to  the  back  during  their  many  trips.  Their  results  were  fitted  to  an 
empirical  equation  to  provide  a  means  of  estimating  food  and  rest  re¬ 
quirements  in  advance  of  any  given  trip.  The  equation  agrees  very  well 
with  the  Data  Book  results  given  above: 

E  -  100(10  +  R  +  2C  +  4H) 

where 

E  =  energy  expenditure  (kcal/day) 

R  =  distance  on  roads  (mi/day) 

C  =  distance  cross-country  (mi/day) 

H  =  height  climbed  (kft/day) 

Energy  expenditure  for  road  travel  is  1/20  kcal/m,  and  energy-rate  expen 
diture  for  a  10  percent  grade  is  twice  that  on  a  horizontal  path  cross¬ 
country. 

Reference  3  indicates  that  4000  kcal/day  is  about  the  maximum  rate 


And  for  other  reasons  discussed  in  the  next  section. 
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of  expenditure  to  be  expected  on  a  continuous  basis  by  a  young,  140-lb 
man  in  good  condition.  If  the  man  is  heavier  than  140  lb  or  is  carrying 
a  backpack,  the  value  of  E  must  be  increased  proportionately. 

One  of  the  first  cost  functions  used  in  this  model  was  of  the  form 
of  Eq,  (1): 


Cost  =  C[1  +  a(grade)  4-  b(grade)^]  (1) 

where  C  -  250  kcal,  the  energy  required  to  travel  between  adjacent  grid 
points  (2.5  km)  on  flat  terrain  (1/10  kcal/m) .  By  choosing  a  4  10b  =  10, 
we  can  double  the  energy  required  for  10  percent  upgrades,  and  by  making 
b  large  or  small  we  can  adjust  the  cost  of  steeper  grades. 

This  form  of  the  cost  function  is  not  as  useful  as  it  appears. 

Since  the  cost  for  the  total  route  is  the  sum  of  the  costs  between  grid 
points  and  grades  are  proportional  to  altitude  changes,  between  the  same 
start  and  end  points  the  linear  term  contributes  the  same  cost  to  all 
routes,  because  the  total  cost  due  to  this  term  depends  only  on  the  alti¬ 
tude  difference  between  the  start  and  end  points.  Similarly,  the  con¬ 
stant  C  will  not  change  the  choice  of  routes. 

A  more  useful  form  of  the  cost  equation  is 

Cost  =  14-  (a|grade|)P  (2) 

By  taking  the  absolute  value  of  the  grade  we  do  not  obtain  the  same  re¬ 
sult  for  different  routes  when  the  power  p  is  1.  Values  of  1  and  2  were 
investigated  for  p;  and  values  of  10,  20,  and  40  for  a.  The  discussion 
of  energy  requirements  suggests  a  value  for  a  =  10  to  double  the  rate 
of  energy  consumption  on  a  10  percent  grade.  This  value  of  a  is  too  low 
because  of  the  interpolation  process. 

Figure  4  is  a  graph  of  the  mean-square  difference  in  altitude  at 
two  grid  points  separated  by  a  distance  £.  Two  curves  are  calculated, 
one  using  interpolated  altitudes  and  another  using  altitudes  read  di¬ 
rectly  from  the  map  at  2.5-km  intervals.  Since  our  interpolation  is  a 
smoothing  process,  we  would  expect  interpolated  altitude  differences  to 
be  less  than  the  actual  differences.  At  the  2.5-km  grid  distance,  the 
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mean- square  difference  is  a  factor  of  5  less  than  the  map  data.  This 
can  be  accounted  for  by  a  factor  of  2  or  more  in  a. 

Furthermore,  since  the  real  terrain  between  grid  points  is  not  a 
smooth  ramp,  extrapolation  of  the  map  data  structure  function  to  smaller 
intervals  indicates  that  the  true  root-mean-square  altitude  difference 
does  not  drop  as  fast  as  the  interval  distance,  and  hence  the  average 
grade  (absolute  value)  is  larger  still.  Using  a  1-km  grid  distance 
would  have  required  increasing  a  by  another  factor  of  2.  We  feel,  there¬ 
fore,  that  a  =  40  is  a  more  realistic  choice  if  the  cost  function  is  to 
be  proportional  to  energy  consumption. 

Because  we  found  so  much  evidence  in  the  literature  and  in  conversa¬ 
tion  with  people  acquainted  with  infiltration  tactics  that  ridge  routes 
and  valley  routes  are  to  be  preferred  to  travel  along  hillsides  or  moun¬ 
tainsides,  we  chose  to  add  a  term  to  the  cost  equation  which  penalized 
grid  points  that  had  steep  gradients,  independent  of  the  relation  of  the 
gradient  direction  to  the  direction  of  travel.  The  term  added  we  call 
the  roughness  factor,  R,  where 

2 

R  =  (3  grade) 

and  the  square  of  the  grade  at  a  grid  point  P  is  found  by  adding  the 
square  of  the  altitude  difference  of  the  two  grid  points  north  and  south 
of  P  to  the  square  of  the  altitude  difference  at  the  two  points  east  and 
west  of  P  and  dividing  by  4  times  the  square  of  the  grid  distance. 

Values  of  3  of  20  and  40  were  used  in  order  that  the  calculated  grade 
using  interpolated  data  would  compare  to  the  real  grade.  This  roughness 
factor  compensates  in  part  for  the  fact  that  the  interpolation  technique 
underestimates  valleys  and  ridges,  as  was  discussed  in  an  earlier  sec¬ 
tion.  Its  use  in  the  cost  equation  is  purely  ad  hoc,  however,  since  few 
data  have  been  found  which  relate  energy  saved  using  ridge  and  valley 
routes  rather  than  hillsides  except  for  the  fact  that  hillsides  will 
often  be  rougher  ground  due  to  erosion.  There  are  important  exceptions, 
however.  Ridgelines  on  new  mountain  ranges  can  be  extremely  rugged. 
Nevertheless,  this  does  not  appear  to  be  the  case  in  Southeast  Asia,v  u 
where  the  mountains  are  older  and  more  eroded. 
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IV,  THE  COMPUTER  MODEL 

PREMIR  is  primarily  a  mathematical  computer  model  which  in  its 
present  form  takes  into  account  terrain  data  from  a  geographical  area 
and  by  means  of  a  suitable  weighting  function  arrives  at  a  set  of  pre¬ 
dicted  (optimal)  infiltration  routes  between  various  start  and  specific 
end  points.  A  total  cost  measure  is  assigned  to  each  route  according 
to  its  computed  degree  of  difficulty.  This  section  will  describe  the 
details  of  the  computer  model  and  how  optimal  routes  are  determined. 

A  basic  underlying  assumption  in  the  development  of  the  model  is 
that  topography,  as  well  as  other  terrain  factors,  is  an  important  de¬ 
terminant  in  the  choice  of  routes.  Thus,  input  to  the  model  includes 
altitude  readings  from  a  chosen  geographical  area  arrayed  in  a  rectan¬ 
gular  grid  format.  A  two-dimensional  interpolation  function  is  employed 
to  interpolate  for  points  between  the  input  grid  data  points,  resulting 
in  a  finer  grid  representation  of  the  area  under  study. 

The  above  procedure  facilitates  the  construction  of  a  somewhat 
comprehensive  network  representation  of  the  study  area  consisting  of 
nodes  and  arcs.  The  nodes  (cross-points)  of  the  network  are  the  input 
and  interpolated  data  points.  Each  interior  point  in  the  network  has 
eight  surrounding  points,  and  each  exterior  point  has  either  three  or 
five  surrounding  points,  depending  on  its  position.  An  example  of  a 
grid  network  is  shown  below;  i  is  an  exterior  point,  and  k,  an  interior 
point. 

o  o  o 

oo.o 

k 

o  o  j  o 

Movement  within  the  network  is  permitted  between  each  point  and  its 
surrounding  points  via  arcs.  Arcs  are  characterized  by  their  beginning 
and  ending  nodes  and  an  associated  cost. 

A  weighting  function  is  used  to  compute  the  cost  of  moving  from 
one  point  within  the  network  system  to  another  as  a  function  of  the 
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change  in  altitude.  Also  available  is  an  option  to  include  in  the 
weighting  function  a  measure  for  the  roughness  of  the  area  within 
which  a  data  point  lies.  All  of  the  data-point  values  (input  and  in¬ 
terpolated)  and  cost  calculations  may  be  printed  as  output  from  the 
model. 

A  dynamic-programming  minimum-path  algorithm  is  then  applied 
to  find  minimum-cost  (optimal)  paths  between  all  possible  start  and 
specific  end  points  of  the  system.  The  algorithm  encompasses  a  value- 
iterative  scheme,  and  the  solution  which  is  produced  gives  the  best 
node  to  go  to  from  each  nodal  point  in  the  system  en  route  to  a  spe¬ 
cific  terminal  point.  Also  obtained  is  a  cost  measure  at  each  nodal 
point  representing  the  total  cost  of  an  optimal  route  from  that  point 
to  the  terminal. 

The  model  is  written  in  FORTRAN  IV  and  under  its  present  config¬ 
uration  allows  for  multiple  cases  to  be  included  during  a  single  com¬ 
puter  run,  subject  to  the  time  constraints  imposed  by  the  operating 
system.  Computer  storage  requirements  are  mainly  dependent  upon  the 
total  number  of  grid  data  points  (both  input  and  interpolated)  for 

the  study  area.  Currently,  a  case  resulting  in  a  100  x  100  grid  ar- 
* 

ray  requires  444k  bytes  of  storage  (k  =  1024)  on  the  IBM  360  system. 
Smaller-sized  problems  require  less  storage  space.  Copies  of  the  com¬ 
puter  source  deck  or  object  deck  (compiled  on  the  optimal  compiler)  may 
be  obtained  from  the  authors. 

PREMIR  is  structured  into  modules  with  each  one  performing  dif¬ 
ferent  tasks.  The  following  is  a  functional  description  of  the  prin¬ 
cipal  modules  (subroutines)  and  underlying  rationale."' 

MAIN  PROGRAM 

The  main  program  functions  as  executor  for  the  model.  Control 
information  relating  to  the  problem  case  as  well  as  the  input  data  on 
the  study  area  is  read  into  the  model  by  this  routine.  It  also  calls 
on  other  routines  to  perform  certain  tasks  as  specified  by  the  user. 

The  main  program  provides  for  special  efforts  to  detect  input  errors 
and  identify  them  for  the  user. 

* 

Such  an  array  contains  10,000  data  points. 
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INTERPOLATION 

The  user  may  choose  to  use  only  the  input  data  or  have  the  model 
interpolate  for  a  variable  number  of  additional  data  points  in  creating 
the  network  representation  of  the  study  area.  In  the  latter  case,  a 
two-dimensional  interpolation  algorithm  is  employed.  The  algorithm 
interpolates  purely  on  the  basis  of  an  inverse-squared  distance  weight¬ 
ing  function.  A  function  of  this  type  not  only  gives  seemingly  satis¬ 
factory  empirical  results,  but  also  presents  an  easy  calculation.  In 
this  function,  the  value  at  any  point,  say  P,  is  a  weighted  average  of 
the  values  of  the  input  data  points  within  a  neighborhood  of  P. 

Here  we  let  Z ^  be  the  value  at  input  data  point  D_^,  and  R^  be  the 
cartesian  distance  between  P  and  D.  The  value  at  point  P,  considering 
N  nearby  input  data  points  somewhat  symmetrically  spaced  about  P,  is 
given  by  the  following: 


f(P)  - 


N  o  N 

^(R.)  Z./^(R.)  for  R^  ^  0  and  all  within 

i  i  the  neighborhood 


Z. 

x 


for  R.  =  0  and  some  D.  within 

1  i 

the  neighborhood 


(3) 


COST  CALCULATION 

In  order  to  obtain  a  set  of  predicted  (optimal)  routes,  a  method 
of  categorizing  all  possible  paths  according  to  their  relative  degree 
of  difficulty  must  be  used.  We  achieve  this  categorization  by  assign¬ 
ing  a  cost  to  travel  between  the  different  nodes  (data  points)  in  the 
network  system  representing  the  study  area.  Thus,  each  allowable  arc 
in  the  system  has  an  associated  cost,  and  the  total  cost  assigned  to  a 
route  from  various  points  in  the  network  to  a  specific  end  point  is 
merely  the  sum  of  the  costs  for  the  arcs  included  in  the  route. 

The  structure  of  PREMIR  facilitates  the  testing  of  various  cost 
functions  which  take  into  account  topography  and  other  relevant  data. 
In  the  present  application,  several  different  cost  function  formats 
and  parameters  were  investigated.  The  model  as  presently  programmed 
employs  a  cost  function  based  on  the  change  in  altitude  between  nodal 
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points  and  an  option  is  available  to  include  a  weighting  for  the  rough¬ 
ness  of  the  area  within  which  a  data  point  lies. 

The  cost  for  traveling  from  one  node,  say  node  i,  within  the  net¬ 
work  system  to  some  adjacent  node  j  is  given  by  the  following: 


Cost  (i,  j) 


1  +  c 


(° 


altitude^  for  two  nodes  which  are 


}  nondiagonally  separated 


1.4  + 


ca 


altitude^ 


for  two  nodes  which  are 
diagonally  separated 


(4) 


where  c,  a,  and  p  are  different  parameters  of  the  cost  function;  G  is 
the  nondiagonal  (horizontal  or  vertical)  distance  in  meters  between 
nodes  of  the  rectangular  network. 

If  the  option  to  include  the  roughness  factor  into  the  cost  calcu¬ 
lation  is  elected,  the  following  procedure  is  used.  Let  (x,  y)  be  the 
cartesian  coordinates  of  node  j.  Then  to  Eq.  (4)  above  we  add  the  fol¬ 
lowing  quantity: 


[^(altitude  (x  +  1,  y)  -  altitude  (x  -  1,  y))/(2G)]2 

2 

+  [B(altitude  (x,  y  +  1)  -  altitude  (x,  y  -  1))/(2G)J  (5) 


where  8  is  a  parameter  of  this  segment  of  the  cost  function. 

To  insure  that  predicted  routes  do  not  extend  beyond  the  grid  ar¬ 
ray  or  go  into  undesired  areas  (e.g.,  into  the  ocean)  a  blocking  tech¬ 
nique  is  used.  That  is,  a  special  altitude  value  is  assigned  to  the 
outside  borders  of  the  grid  by  another  subroutine  of  the  model.  If  a 
section  of  the  rectangular  grid  includes  a  large  body  of  water  or  other 
undesired  area  of  travel,  the  altitude  values  (data  points)  in  the  sec- 

tion(s)  can  also  be  flagged.  In  computing  the  costs,  a  check  is  made, 

6 

and  if  the  special  altitudes  are  sensed,  a  very  high  cost  (10  )  is  as¬ 
signed  for  traveling  to  that  point.  An  input  to  the  model  is  an  upper 
bound  for  the  special  altitudes,  and  all  points  with  altitudes  less 
than  or  equal  to  this  bound  are  assigned  the  high  cost. 

What  effect  blockage  of  certain  passageways  will  have  on  a  set  of 
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predicted  routes  may  be  investigated  by  inputting  to  the  model  coordi¬ 
nates  of  the  place(s)  to  be  blocked  and  assigning  a  high  cost  accord- 
ingly.  However,  the  cost  calculations  for  neighboring  points  are  also 
affected  when  blocking  is  employed. 

DYNAMIC  PROGRAMMING 

A  network  structure  consisting  of,  say,  p  nodes  with  directed 
arcs  connecting  some  of  the  nodes  and  costs  assigned  for  traversing  the 
various  arcs  is  amenable  to  solution  by  the  techniques  of  dynamic  pro¬ 
gramming.  The  term  dynamic  programming  is  frequently  applied  to  the 
mathematical  analysis  of  problems  in  which  conditions  that  must  be  sat¬ 
isfied  by  an  optimal  time-staged  decision  process  are  to  be  exploited 
to  determine  best  action. 

PREMIR  incorporates  a  dynamic -programming  minimum- cost -path  algo¬ 
rithm  to  find  optimal  routes  through  the  network.  The  solution — an 
optimal  solution--is  defined  to  be  a  path  (from  nodal  points  in  the 
network  to  a  specific  terminal  point)  that  satisfies  the  following  prin¬ 
ciple  of  optimality: 

Regardless  of  the  previous  rationale  used  to  arrive 
at  a  particular  state ^  the  remaining  decisions  as 
to  what  path  to  take  en  route  to  the  terminal  must 
themselves  constitute  an  optimal  solution . 

Here  we  let  be  the  present  value  (sum  of  the  costs)  of  an  op¬ 
timal  path  from  node  i  to  the  terminal  node  r,  where  optimal  implies  a 
path  having  minimal  total  cost.  Now,  if  an  optimal  path  from  node  i 
to  node  r  starts  by  first  going  to  node  j,  then 


and 


+  C.  . 
iJ 


Y1  +  CM 

k  xk 


(6) 

(7) 


for  all  k  4  j  in  the  system. 
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Equation  (6)  states  that  for  an  optimal  path  from  node  i  to  the 
terminal,  r,  which  starts  by  going  to  node  j,  the  cost  is  the  sum  of 
the  cost  of  an  optimal  path  from  node  j  to  r  and  the  cost  of  going  from 
node  i  to  node  j.  Equation  (7)  insures  that  among  the  possible  choices 
of  paths  to  take  from  node  i,  none  is  better  than  the  one  starting  with 
node  j  . 

Since  Eqs.  (6)  and  (7)  must  hold  for  every  node  i  different  from 
the  terminal  node  r,  then  must  satisfy  a  set  of  functional  equations, 

Y ,  =  minimum  [ Y .  +  C ,  .  ] 

1  «,  j>  j 

in  the  system  (8) 

for  all  i  4  r,  and 

Y  e  0  (9) 

r 

A  method  of  successive  approximation  may  be  used  to  solve  for  the 
p  -  1  unknowns,  Y_^(i  4  r)  .  Specifically,  a  value-iterative  algorithm 
is  applied.  The  value  of  the  Y^.!s  is  initially  set  at  zero,  and  the 
Y^'s  are  then  calculated  as  in  Eq.  (8)  for  every  node  i  4  r  in  the  sys¬ 
tem.  That  is,  the  first  iteration  is  computed  as  follows: 

Y (1,  i)  =  minimum  [0  +  C.  .] 

(i,  j)  1J 

in  the  system  (10) 


for  all  i  4  r  with  Y(l,  r)  =  0.  The  first  index  on  Y  is  the  iteration 
number. 

th 

On  the  n  iteration,  the  Y^! s  computed  during  iteration  n  -  1 
become  estimates  of  the  Y^'s,  and  new  Y^*s  are  computed  as  follows: 


Y(n,  i)  =  minimum  [Y(n  -  1,  j)  +  C .  .  ] 

(i,  j)  J 

in  the  system 


(11) 


again  for  all  i  4  r  and  Y(n,  r)  =  0. 
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The  quantity  Y(n,  i)  is  interpreted  as  the  minimum  cost  of  a  path 

starting  from  node  i  that  contains  exactly  n  arcs,  unless  the  terminal 

is  on  the  path,  in  which  case  the  path  terminates  at  node  r.  If  the 

k 

cost,  C..,  around  every  loop  in  the  system  is  positive,  then  an  n 
exists  such  that  for  all  n  ^  n  ,  Y(n,  i)  =  Y (n  ,  i) .  We  terminate  the 
iterative  process  when  this  condition  is  met.  That  is,  after  each  it¬ 
eration  greater  than  some  q(q  >  1) ,  we  test  Y(n,  i)  =  Y(n  -  1,  i),  and 
when  this  condition  is  met  for  all  ifs  in  the  system,  we  have  obtained 
the  optimal  solution.  We  also  obtain  another  vector,  g(n,  i) ,  which 
has  its  value  set  equal  to  the  j  that  minimizes  Eq.  (10)  or  (11),  de¬ 
pending  on  the  iteration  number. 

That  the  process  will  terminate  in  a  finite  number  of  steps  is 
shown  by  the  following  argument:  As  n  increases,  eventually  every 
path  of  n  arcs  reaches  the  terminal.  And  when  this  occurs,  Y(n,  i) 
will  stabilize  to  the  correct  minimal-cost  value  for  the  path,  since 
the  only  new  value  added  to  Y(n,  i)  will  be  Y(n,  r)  =0. 

i 

There  are  several  algorithms  for  finding  optimal  paths  through 
a  network.  Corresponding  to  almost  any  shor test-path  algorithm  some 
class  of  network  structures  exists  for  which  the  algorithm  is  relatively  1 
inefficient.  Since  all  iterative  procedures  are  in  some  sense  ineffi¬ 
cient,  economic  considerations  alone  demand  that  care  be  taken  in  the 

selection  and  adaptation  of  a  shortest-path  algorithm. 

(4) 

Dreyfus  has  examined  various  shortest-path  algorithms  by  devel¬ 
oping  upper  bounds  on  the  maximum  number  of  computational  steps — the 
number  of  additions  and  comparisons  directly  attributable  to  the  algo¬ 
rithm  itself — involved  in  solving  a  standardized  network  problem.  He 
considers  one  procedure  to  be  significantly  superior  to  another  when 
the  bounds  differ  by  a  multiplicative  factor  involving  p,  the  number 
of  nodes.  By  this  criterion  the  algorithm  applied  in  the  present  study 
is  less  efficient  than  some  other  applicable  algorithms.  That  is,  the 
algorithm  described  above  has  a  higher  bound  on  the  maximum  number  of 
computational  steps  than  some  other  algorithm  for  solving  the  standard¬ 
ized  problem.  However,  this  assertion  and  its  implications  can  be 

Such  algorithms  are  referred  to  as  shortest-path  algorithms  where 
the  length  of  each  (i,  j)  is  taken  to  be  the  cost  . 
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mis  leading,  since  the  actual  number  of  computational  steps  required  to 
find  solutions  in  this  investigation  falls  far  short  of  the  bound  com¬ 
puted  by  Dreyfus  and  the  actual  number  of  steps  required  by  some  other 
algorithms.  (This  result  is  in  part  due  to  certain  specifics  of  the 
problem  and  modifications  of  the  algorithm  for  programming  purposes.) 

From  our  limited  experience  in  this  investigation  it  appears  that 
the  actual  number  of  computational  steps  required  is  more  dependent 
upon  the  interrelation  among  the  costs,  C..,  than  upon  the  upper  bounds. 
Also  dynamic-programming  algorithms  can  require  a  great  deal  of  com¬ 
puter  storage  and  time  to  solve  various  problems.  Sometimes  the  size 
requirements  of  a  particular  problem  may  be  such  that  the  modifications 
necessary  to  adapt  a  specific  algorithm  may  themselves  render  the  algo¬ 
rithm  inefficient.  The  inability  to  specify  the  actual  number  of  com¬ 
putational  steps  required  is  a  shortcoming  of  many  shortes t-path 
algorithms. 

In  real  applications  of  PREMIR,  many  initial  and  terminal  points 
will  be  considered  and  it  may  well  be  inefficient  to  make  a  separate 
run  for  each  case.  Also,  in  a  more  rigorous  study,  different  varia¬ 
tions  of  the  problem  posed  here  may  be  investigated,  such  as  determin- 
th 

ing  the  K  best  route  or  finding  the  best  route  visiting  specified 
nodes.  Dreyfus  also  describes  algorithms  developed  by  Floyd  and 
Dantzig  for  finding  simultaneously  the  shortest  paths  between  all 
pairs  of  nodes  of  a  network,  algorithms  for  finding  the  second-  (and 
third-,  etc.)  shortest  path  through  a  network,  and  a  method  for  find¬ 
ing  shortest  paths  visiting  specified  nodes.  Given  the  necessary  time 
and  resources,  worthwhile  experience  with  the  model  could  be  gained  by 
exploration  with  some  of  these  ideas  as  options  or  main  features. 


Is' 
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V*  THE  VIETNAM  EXAMPLE 


Since  we  are  limited  by  computer  size  to  approximately  10,000  grid 
points,  we  can  choose  to  examine  a  large  area  with  a  coarse  grid  or  a 
smaller  area  with  a  fine  grid.  The  proper  order  would  seem  to  be  to 
find  the  optimal  routes  over  the  larger  area  and  examine  details  of  the 
routes  at  a  later  time  with  a  finer  grid. 

The  area  chosen  is  illustrated  in  Fig.  5.  The  coordinates  closest 
to  the  map  are  the  UTM  coordinates  of  the  10-km  grid  network.  The  sec¬ 
ond  row  of  coordinates  (at  the  very  bottom  and  far  left)  are  computer 
coordinates  established  for  this  problem.  One  unit  in  these  coordinates 
represents  a  distance  of  2.5  km,  the  grid  size  used  to  establish  routes. 
We  shall  use  these  coordinates  in  the  discussion  which  follows. 

The  computer  run  requires  that  a  terminal  point  be  established; 
minimum-cost  paths  from  every  other  point  on  the  grid  are  then  computed. 
Since  the  optimum  path  from  any  point  is  unique  and  independent  of  the 
path  to  that  point,  it  is  sufficient  for  the  computer  to  print  out  a 
table  which  specifies  to  which  of  the  eight  adjacent  points  the  user 
proceeds  from  any  point  on  the  grid.  From  this  table  optimal  paths  to 
the  terminal  point  can  be  traced  from  any  initial  point. 

In  Fig.  5  the  terminal  point  chosen  was  at  coordinates  (60,  66), 
the  east/west  coordinate  being  given  first.  Initial  points  were  chosen 
every  10  km  along  the  left-hand  edge  of  the  grid.  The  cost  function  of 
Eq.  (1)  was  used  with  a  =  10  and  b  =  50.  Since  this  set  of  values  under¬ 
estimates  the  costs  of  rough  terrain,  we  see  very  little  route  selection 
to  avoid  rough  terrain. 

In  Fig.  6  we  have  increased  the  cost  parameters  to  a  =  20  and 
b  =  200.  We  can  see  now  some  further  terrain  avoidance,  but  the  routes 
still  tend  to  be  rather  straight. 

Since,  as  discussed  before,  the  linear  term  of  cost  Eq.  (1)  pro¬ 
vides  no  discrimination  on  route  choice,  we  performed  the  remainder  of 
our  calculations  with  cost  Eq.  (2).  Figure  7  illustrates  such  a  case 
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where  the  parameters  have  their  smallest  value,  p  =  1  and  a  =  10.  The 
terminal  point  is  at  (64,  47),  a  likely  base  area,  and  the  initial 
points  were  chosen  to  lie  near  to  the  road  network  in  Laos.  The  rea¬ 
son  for  these  initial  points  is  that  the  infiltrators  can  in  theory 
be  brought  to  these  points  by  truck,  and  the  costs  of  reaching  these 
points  would  be  on  a  different  basis  than  the  costs  for  reaching  them 
on  foot. 

In  Fig.  7  there  is  little  terrain  discrimination  in  the  route  se¬ 
lection  because  the  value  of  a  underestimates  the  cost  of  rough  terrain. 
In  Fig.  8,  however,  we  have  increased  a  by  a  factor  of  4  to  its  real¬ 
istic  value  and  we  have  added  the  roughness  term  to  the  cost  equation 
with  a  value  6  =  40.  Furthermore,  grades  larger  than  10  percent  are 
more  heavily  discriminated  against  by  setting  p  =  2  so  that  terrain 
costs  increase  as  the  square  of  the  grade. 

In  this  figure  we  see  considerable  selectivity  in  routes.  The 
northern  routes  reach  into  Vietnam  to  pass  along  the  A  Shau  valley. 

The  central  route  proceeds  in  a  way  that  is  quite  surprising  in  view 
of  the  contour  lines  shown  on  the  map.  Especially  strange  is  its  cross¬ 
ing  the  valley  at  coordinates  (44,  44).  An  examination  of  the  terrain 
data  and  their  interpolation  indicates  that  our  input  grid  was  too 
coarse  for  this  narrow  valley  and  it  does  not  appear  in  the  interpola¬ 
tion. 

In  Figs.  9  and  10  we  examine  another  possible  terminal  point  at 
(70,  18)  just  inside  the  Laos  border.  In  these  examples  we  can  see  the 
changes  produced  by  adding  the  roughness  factor  (8  =  40)  to  the  cost 
equation  for  Fig.  10.  In  both  figures,  a  =  40  and  p  =  2.  The  effect 
of  the  factor  is  to  shift  many  of  the  central  routes  to  that  path  which 
passes  along  the  Vietnam  border.  The  remaining  central  path  shifts  to 
cross  the  rough  terrain  at  a  narrower  point.  The  southerly  routes  also 
pick  their  way  more  carefully  through  the  rough  terrain.  Although  there 
is  flat  terrain  available  inside  Vietnam,  the  distances  involved  are 
such  that  the  routes  remain  for  the  most  part  inside  Laos  in  spite  of 
the  mountainous  areas . 


Southwest  of  grid  coordinates  (58,  58). 


In  Figs.  11  and  12  we  examine  routes  to  Hue  (72,  67)  with  and  with¬ 
out  the  roughness  factor.  On  the  basis  of  altitude  costs  alone  there  is 
considerable  choice  of  alternate  routes  through  the  mountains,  as  shown 
in  Fig.  11.  Adding  the  roughness  cost  in  Fig.  12  narrows  the  choice  con 
siderably.  In  both  cases  the  central  routes  tend  to  pass  through  the 
A  Shau  valley  in  the  region  of  (56,  56).  Many  of  the  routes  in  Fig.  12 
move . considerably  to  the  north  to  avoid  the  mountainous  regions  in  the 
center . 

This  effect  is  also  present  in  Fig.  13  where  the  terminal  point  has 
been  moved  further  west  to  (60,  66).  The  general  route  structure  is  un¬ 
changed. 

Because  of  the  isolation  of  the  central  route  in  Fig.  13,  it  was 
decided  to  investigate  the  effect  of  a  barrier  built  across  this  route. 
To  build  a  barrier  in  the  model  it  is  sufficient  to  raise  the  costs  at 
a  continuous  set  of  grid  points.  This  was  done  in  Fig.  14  along  a  di¬ 
agonal  line  shown  passing  through  (46,  46).  The  effect  is  to  cause  the 
central  route  to  move  considerably  to  the  north  of  the  original  route 
rather  than  simply  move  around  the  edge  of  the  barrier.  In  this  way 
the  effectiveness  of  various  barrier  locations  and  the  minimum  size  of 
a  barrier  necessary  to  require  a  major  route  shift  can  be  investigated. 

In  Fig.  15  we  see  the  same  case  but  with  cross-shaped  barriers 
placed  in  the  Khe  Sanh  and  Quang  Tri  areas.  Because  of  the  flatness  of 
the  terrain  in  these  areas,  the  infiltration  routes  bypass  these  bar¬ 
riers  with  little  change. 


Fig- 1 1 — Predicted  infiltration  routes,  with  terminal  point  at  coordinates  (72  ,  67)  v  - ... 

(p  ~  2,  a  -  40,  (3  =  0) 


Fig .13 — Predicted  infiltration  routes,  with  terminal  point  at  coordinates  (60,  66) 

(p  =  2,  a  =  40,  £  =  40) 
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Fig.  14 — Predicted  infiltration  routes,  with  diagonal  barrier 
and  terminal  point  at  coordinates  (60,  66) 

(p  -2,  a  =  40,  (3  =  40) 
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VI,  APPLICATIONS  OF  THE  MODEL 

INFILTRATION  ROUTE  CHARACTERISTICS 

It  has  never  been  seriously  suggested  that  infiltration  routes  are 
chosen  at  random.  Traditional  communication  routes  almost  always  fol¬ 
low  easy  rather  than  difficult  terrain.  Even  if  a  route  is  discovered 
accidentally  and  further  exploration  stops,  it  will  have  some  "shortest- 
pathn  characteristics  even  though  the  route  may  not  be  optimum.  The 
model  provides  a  systematic  approach  for  investigating  the  criteria  and 
degree  of  randomness  involved  in  route  determination,  provided  some  in¬ 
filtration  routes  are  known  to  the  investigator. 

For  example,  if  a  given  route  appears  to  follow  easy  terrain,  it 
is  possible  to  use  the  model  to  predict  various  routes  between  end 
points  of  the  known  route,  using  different  kinds  of  terrain  cost  func¬ 
tions.  If  the  real  route  departs  systematically  from  the  "terrain  opti¬ 
mal"  routes  in  some  region,  an  attempt  would  then  be  made  to  determine 
how  the  real  and  predicted  routes  in  this  region  differ.  Questions 
might  be  asked  such  as,  Does  the  real  route  avoid  villages  that  the 
predicted  route  approaches?  Does  the  real  route  detour  to  follow  a  tra¬ 
ditionally  established  path?  Does  the  real  route  follow  a  path  where 
vegetation  provides  cover  from  the  air?  Could  the  deviation  be  due  to 
the  existence  of  an  unknown  base  camp  or  logistics  camp?  Once  a  reason 
for  the  deviation  is  conjectured,  a  new  cost  function  may  be  constructed 
which  includes  a  weighting  factor  for  the  new  condition,  and  the  process 
is  repeated  using  new  route  predictions  and  looking  for  further  depar¬ 
tures  from  the  predicted  routes. 

It  is  clear  from  this  discussion  that  route  characteristics  other 
than  the  optimizing  criteria  can  be  inferred.  For  example,  it  may  be 
possible  to  discover  camps  or  rallying  points  by  examining  the  deviation 
of  real  routes  from  those  predicted.  When  fragmentary  route  information 
exists,  it  may  be  possible  to  determine  start  or  end  points  for  the  real 
routes  by  matching  the  route  fragments  to  routes  predicted  using  an  ar¬ 
ray  of  start  and  end  points. 

By  use  of  sensors  along  known  and  predicted  routes,  it  may  be  pos¬ 
sible  to  determine  shifts  of  traffic  from  one  route  to  another;  and  by 
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use  of  the  model,  to  determine  whether  the  traffic  shift  is  due  to  a 
change  of  objective  (end  point)  or  a  change  of  selection  criterion  (for 
example,  heavy  logistics  require  easier  terrain) . 

The  alternate-route  spectrum  can  also  be  investigated  by  the  model. 
For  each  predicted  route  we  can  examine  alternate  routes  which  have  a 
cost  only  slightly  larger  than  that  of  the  optimum  route.  In  the  case 
where  small  changes  in  the  route  result  in  large  increases  in  route  costs 
we  can  be  confident  the  predicted  route  is  significant.  When  many  alter¬ 
nate  routes  exist  at  only  slight  cost  increases,  then  the  predicted 
route  should  not  be  taken  too  much  for  granted  and  the  whole  alternate- 
route  array  should  be  examined. 

BARRIER  AND  SENSOR  EMPLACEMENT 

In  many  cases  it  may  be  necessary  to  construct  barriers  or  strong 
points  to  prevent  or  restrict  infiltration.  Since  the  resources  avail¬ 
able  to  the  infiltrator  are  usually  limited,  some  kind  of  criterion  for 
allocating  resources  for  barrier  construction  is  desirable.  If  the  ter¬ 
rain  is  at  all  rugged,  it  would  not  seem  wise  to  construct  a  barrier 
using  the  criterion  of  equal  cost  per  mile  of  border.  More  resources 
should  be  expended  where  the  cost  to  the  infiltrator  is  otherwise  small. 

One  approach  to  barrier  resource  allocation,  for  example,  would  be 
to  prepare  a  route  map  from  the  model,  such  as  that  shown  in  Fig.  13, 
using  a  number  of  start  points  and  one  or  more  appropriate  end  points. 
Where  several  routes  join  to  pass  through  choke  points,  barriers  can  be 
emplaced  in  the  model  as  was  done  in  Fig.  14,  and  the  increased  cost  of 
the  new  route  can  be  examined.  If  the  increase  is  sufficiently  large, 
a  fairly  substantial  barrier  would  be  indicated. 

A  more  systematic  approach  would  be  to  examine  infiltration  costs 
along  the  entire  border  (or  barrier  line)  by  choosing  points  along  the 
border,  first  as  terminal  points  from  and  then  as  start  points  for  some 
point  within  the  defended  country,  and  adding  the  costs  at  each  border 
point  for  the  two  routes  leading  to  and  away  from  the  point.  Since  the 
barrier  represents  an  additional  cost  to  the  infiltrator  and  this  bar¬ 
rier  penetration  cost  should  increase  as  the  cost  of  the  barrier  in¬ 
creases,  then  barrier  resources  can  be  allocated  so  that  the  total  cost 
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to  the  infiltrator  (infiltration  plus  barrier  penetration)  is  the  same 
for  all  points  of  border  crossing.  .J 

An  extension  of  this  approach  allows  the  determination  of  the  op¬ 
timum  location  of  the  barrier  line.  When  the  barrier  along  a  given  line 
is  optimized,  the  infiltration  cost  is  uniquely  determined.  Holding  the 
resources  for  barrier  emplacement  fixed,  and  recognizing  that  barrier 
costs  increase  as  the  barrier  is  placed  further  into  the  mountains  away 
from  populated  areas,  we  can  determine  where  the  barrier  line  for  fixed 
total  resources  presents  the  highest  total  cost  to  the  infiltrator. 

This  approach  can  also  be  used  for  the  placement  of  sensor  fields. 

It  is  assumed  that  the  infiltrator  would  like  to  avoid  sensor  fields; 
increasing  the  density  of  these  fields  in  regions  of  easiest  travel, 
thereby  causing  the  infiltrator  to  avoid  these  regions,  will  increase 
his  costs. 

Identification  of  choke  points  (i.e.,  points  around  which  alternate 
routes  have  a  large  incremental  cost)  suggests  positions  for  strong-point 
defenses.  It  was  evident  without  analysis  why  Khe  Sanh  was  a  good  loca¬ 
tion  for  restricting  infiltration.  Quite  objectively,  Fig.  13  points  up 
distinctly  the  value  of  the  Route  9  corridor.  However,  the  model  also 
indicates  the  width  of  the  corridor  and  the  range  of  routes  over  which 
infiltration  costs  are  relatively  stable.  To  be  effective,  a  strong 
point  at  Khe  Sanh  would  have  to  defend  a  very  large  area. 

OFFENSE  APPLICATIONS 

The  model  should  be  of  value  to  the  infiltrator  as  well  as  the  de¬ 
fender.  It  can  recommend  routes  quickly  in  areas  where  the  infiltrating 
force  has  little  experience.  Perhaps  the  most  effective  criticism  of 
the  model  used  by  the  defender  is  that  the  infiltrator  may  not  be  as 
systematic  as  the  model  implies.  This  problem  is  not  important  to  the 
infiltrator  using  the  model.  He  can  estimate  the  costs  of  terrain  for 
known  personnel  and  equipment.  He  knows  his  start  and  end  points.  He 
knows  how  to  weight  population  factors  and  detection  by  air,  etc.  Except 
for  having  to  consider  unknown  hazards  and  factors  which  would  be  a  prob¬ 
lem  in  any  case,  the  model  can  produce  a  number  of  preferred  infiltra¬ 
tion  routes  and  estimates  of  the  costs  and  time  required  to  traverse 
each  route. 
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Other  applications  of  the  model  can  be  suggested;  in  many  cases 
these  are  similar  to  techniques  already  being  applied.  Choosing  a 
minimum- cost  route  to  construct  a  road  between  two  points  is  an  obvi¬ 
ous  application.  Finding  the  best  routes  for  ground  resupply  of  an 
advanced  post  is  another.  Finding  the  most  effective  ambush  points 
when  a  sensor  signals  traffic  at  some  distance  is  a  possibility.  There 
exists  a  range  of  applications  such  as  determining  boundaries  of  polit¬ 
ical  subdivisions  on  the  basis  of  ease  of  communication  which  is  too 
broad  for  further  consideration  at  this  time. 

Considerably  greater  experience  with  the  model  is  required  to 
fully  estimate  its  capabilities  and  potentialities.  Its  concept  is 
simple,  but  the  diversity  of  parameters  it  encompasses  provides  a  rich 
field  for  further  investigation. 
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Appendix  A 

A  JOSS  PROGRAM  FOR  INVERSE-SQUARE  INTERPOLATION 

When  data  are  available  on  a  two-dimensional  rectangular  grid,  it 
is  often  of  interest  to  determine  values  and  contour  lines  on  a  scale 
finer  than  the  data  grid  by  interpolation.  There  are  many  ways  inter¬ 
polation  can  be  done,  but  we  chose  an  inverse-square  technique.  We 
present  here  a  simple  JOSS  program  which  can  be  adapted  to  other  remote 
console  systems  with  similar  formatting. 

The  program  listed  below  accepts  49  data  inputs  (from  a  7  x  7  grid) 
as  z(i,  j).  Parts  1,  2,  and  3  input  the  data.  Step  3.2  was  added  to 
avoid  having  to  input  a  decimal  point.  Part  99  can  be  used  to  start 
the  program  if  the  data  are  already  in  storage. 

The  interpolated  value  at  each  point  is  calculated  as  a  weighted 
average  of  the  16  nearest  input  grid  values.  The  weighting  function  is 
proportional  to  the  square  of  the  reciprocal  of  the  distance  between 
the  interpolation  point  and  each  of  the  16  data  points.  The  calcula¬ 
tion  proceeds  by  interpolating  a  horizontal  line  at  the  top  of  the  area, 
printing  that  line,  and  then  repeating  the  process  for  the  next  line 
down,  and  so  on. 

The  interpolated  area  lies  within  the  5x5  subgrid  of  the  origi¬ 
nal  7x7  data  grid.  We  have  chosen  as  an  example  an  interpolation  in¬ 
terval  of  one- tenth  of  the  original  data-grid  interval. 

Step  1.2  does  the  vertical  scan  of  the  program.  Step  4.1  activates 
parts  5,  8,  9,  and  7  to  calculate  the  interpolation  along  a  horizontal 
line.  Each  interpolated  value  is  calculated  as  h(x)  in  Step  5.91. 

Steps  4.2  and  4.3  then  initiate  an  unlikely  looking  set  of  calcu¬ 
lations  in  parts  10  and  11,  which  arise  because  of  the  formatting 
limitations  of  JOSS.  Part  10  selects  one  digit  to  be  printed  at  the 
interpolated  point  representing  a  truncated  value  of  h(x)  .  Part  11 
combines  five  of  these  digits  into  a  single  number  b(k),  alternating 
the  digits  with  the  digit  1,  which  is  a  dummy.  The  purpose  of  the  dummy 
is  to  space  the  digits  horizontally  equal  to  the  vertical  spacing  of 
the  printed  lines.  Since  JOSS  will  not  print  a  leading  zero  in  a  num¬ 
ber,  Step  11.3  was  added  to  use  a  dash  to  avoid  an  empty  space  where  a  /-• 
zero  should  be. 
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Steps  4.31,  4.32,  and  4.33  do  the  same  as  part  11  except  for  com¬ 
bining  four  digits  rather  than  five,  since  we  could  not  fit  the  fifth 
digit  into  the  final  space  of  the  form  line.  Step  4.4  does  the  final 
output  for  each  horizontal  line  by  typing  the  eight  numbers,  b(k),  into 
the  form,  giving  a  continuous  line  of  39  digits  separated  by  38  dummy 
digits . 

Figure  A-l  shows  an  area  of  South  Vietnam  drawn  from  a  map  with 
300-m  contour  lines.  Figure  A- 2  is  a  sample  output  from  the  interpola¬ 
tion  program  of  the  same  area,  where  each  digit  (not  a  dummy)  represents 
the  altitude  in  100-m  intervals.  The  data- input  grid  is  the  same  as 
that  shown  in  Fig,  A-l,  but  grid  points  surrounding  the  interpolated 
area  are  not  shown. 
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Fig  .A—  1 — Map  contours 
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Use  file  270  (35350)* 

Roger. 

Recall  item  25  (MAP). 

Done . 

Type  a 11, size. 

1.1  Do  part  2  for  j=5(-l)-l. 

1.15  Page, 

1.16  Set  c=100. 

1.2  Do  part  4  for  y=39(-l)0. 

2.1  Do  part  3  for  i=- 1(1)5. 

3. 1  Demand  z(i , j  ). 

3.2  Set  z(i , j )=z(i , j  )/10 . 

4.1  Do  part  5  for  x=0(l)39. 

4.2  Do  part  10  for  x=0( 1)39. 

4.3  Do  part  11  for  k=0(5)30. 

4.31  Set  b(35)=a(38)+c* (a( 37 )+c» (a(36 )+c*a(35 ) ) ) . 

4.32  Set  b ( 35 )=b( 35  )+1010 10 . 

4.33  Set  b ( 35  )=-b ( 35 )  if  a(35)  =  0. 

4.4  Type  b(0 ) ,b(5 ) ,b( 10 ) ,b ( 15 ) ,b (20 ) ,b(25) ,b( 30 ) ,b ( 35 )  in  form  1. 

5.1  Do  part  8  if  fp(x/10)=0  and  fp(y/l0)=0. 

5.2  Done  if  fp(x/10)=0  and  fp(y/10)=0. 

5.3  Set  N=0. 

5.4  Set  g(x)=0. 

5.5  Set  X=  ip( x/10  ) . 

5.51  Set  Y=ip (y/10  ) . 

5.9  Do  part  9  for  i=X-l(l)X+2. 

5.91  Set  h(x)=g(x)/N. 

7.1  Set  f ( i , j )=  (x/ 10-i )*2t (y/10-i  )*2 . 

7.2  Set  N=N+ 1/f (i , j  ) . 

7.3  Set  g(x)=g(x )+z(i , j )/f (i , j ) . 

8.10  Set  m=x/10. 

8.11  Set  n=y/10 . 

8.2  Set  h(x)=  z(m>n). 

8.3  Done. 

9.1  Do  part  7  for  j=Y-l(l)Y+2. 

10.1  Set  a(x)=ip(10»fp(h(x)/10) ). 

11.1  Set  b(k )=a(k+4 )+c* (a(k+3 )+c*(  a(k+2)+c* (  a(k+l )+c*a(k  ) ) ) ) . 

11.2  Set  b ( k )  =  b(k)+  10101010. 

11.3  Set  b(k)=-b(k)  if  a(k)=0. 

99  To  step  1.15. 

Form  1: 

___  1111  1  1 


size  =  218 


Str<  f  V~. 


J/  L\  j* 
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Appendix  B 

PREPARATION  OF  AN  INPUT  DECK  FOR  THE  FORTRAN  IV 
VERSION  OF  PREMIR 

INPUT  DATA 


Control  Variables 

1.  NtXMX  is  the  number  of  the  last  column  of  data  in  the  input 

matrix  (grid) . 

2.  NUMY  is  the  number  of  the  last  row  of  data  in  the  input 

matrix  (grid) . 

3.  NUMX2  is  the  actual  number  of  columns  of  data  from  the  in¬ 

put  matrix  which  will  be  retained  in  the  subsequent 
network . 

4.  NUMY 2  is  the  actual  number  of  rows  of  data  from  the  input 

matrix  which  will  be  retained  in  the  subsequent  net¬ 
work  . 

5.  NX  is  the  number  of  the  first  column  of  data  in  the  input 

matrix. 

6.  NY  is  the  number  of  the  first  row  of  data  in  the  input 

matrix. 

7.  ND0X  is  the  beginning  horizontal  (X)  index  for  nodes  in  the 

subsequent  network  system. 

8.  ND0Y  is  the  beginning  vertical  (Y)  index  for  nodes  in  the 

subsequent  network  system. 

9.  INCX  is  the  incremental  distance  between  adjacent  nodes  of 

the  network  which  are  horizontally  separated. 

10.  INCY  is  the  incremental  distance  between  adjacent  nodes  of 

the  network  which  are  vertically  separated. 

11.  IC0ST  indicates  whether  or  not  the  input  and  interpolated 

data  values  and  cost  calculations  are  to  be  printed 
as  output.  If  IC0ST  >  0,  the  information  will  be 
printed. 

12.  IBL0CK  is  the  number  of  nodes  in  the  network  to  which  travel 

from  surrounding  nodes  is  to  be  blocked. 

13.  INTER  indicates  whether  or  not  interpolation  is  to  be  applied. 

If  INTER  >  0,  interpolation  will  not  be  applied. 


14,  NDEC 

15.  IR 


16.  IGRIDX 


17.  IGRIDY 


is  the  number  of  decimal  places  to  which  interpolated 
values  are  to  be  rounded. 

indicates  whether  or  not  the  roughness  factor  is  to  be 
included  in  the  cost  calculations.  Input  IR  >  0  for 
inclusion. 

is  the  number  of  interpolated  data  points  +1  between 
adjacent  input  data  points  measured  in  the  horizontal 
direction. 

is  the  number  of  interpolated  data  points  +1  between 
adjacent  input  data  points  measured  in  the  vertical 
direction. 


Cost  Data 

i- 

1.  ZER0  is  the  altitude-limiting  variable.  All  altitude  values 

(input  and  interpolated)  less  than  or  equal  to  ZER0  will 
be  considered  undesirable  for  travel  (e.g.,  the  ocean), 
and  a  high  cost  will  be  assigned  to  them. 


2.  C,  G,  H,  are  parameters  of  the  cost  function  (see  Section  IV). 
ALPHA,  P,  Input  H  as  "1.0". 

BETA 


Dynamic  Program  Data 

1.  ID0X,  IT0X  are  the  horizontal  component  (X)  indexes  of  the  boun¬ 
dary  nodes  (lower  and  upper,  respectively)  for  the  area 
of  the  network  in  which  optimal  routes  are  to  be  pre¬ 
dicted  . 


2.  JD0Y ,  JT0Y  are  the  vertical  component  (Y)  indexes  of  the  boundary 

nodes  (lower  and  upper,  respectively)  for  the  area  of 
the  network  in  which  optimal  routes  are  to  be  predicted. 

3.  I START  is  the  number  of  iterations  the  model  performs  before 

testing  for  the  optimal  solution. 


4.  T(1),'T(2)  are  the  X  and  Y  coordinates  (respectively)  of  the  node 

designated  as  the  terminal. 


CARD  INPUT  F0RMAT 
A  denotes  alphanumeric  data. 

I  denotes  integer  data  right-adjusted  in  its  field. 


F.d  denotes  real  numeric  data  with  d  =  number  of  decimal  points. 
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First  Card 

Col  1-72  Code  name  for  run  (may  be  left  blank). 

Matrix  and  Network  Data 

If  interpolation  is  applied,  then  the  outside  bordering  input  data 
points  may  not  be  included  as  nodes  in  the  network  system,  since  an  inter¬ 
polated  data  point  requires  that  a  somewhat  symmetric  neighborhood  of 
input  data  points  be  used.  The  input  data  are  considered  as  a  mxn  matrix 
and  the  user  inputs  to  the  model  information  relating  to  the  form  and 
structure  of  the  matrix.  The  resulting  network  is  rectangular  in  shape 
with  each  node  having  x,  y  coordinates. 


Second  Card 

Col 

1-5 

Last  column  number  of  input  data  matrix  (NUMX) 

(i) 

Col 

6-10 

Last  row  number  of  input  data  matrix  (NUMY) 

(i) 

Col 

11-15 

Number  of  data  columns  retained  in  network  (NUMX2) 

(i) 

Col 

16-20 

Number  of  data  rows  retained  in  network  (NUMY2) 

(i) 

Col 

21-25 

First  column  number  (NX) 

(i) 

Col 

26-30 

First  row  number  (NY) 

(i) 

Col 

31-35 

Beginning  X  index  for  network  (ND0X) 

(i) 

Col 

36-40 

Beginning  Y  index  for  network  (ND0Y) 

CD 

Col 

41-45 

Horizontal  incremental  distance  (INCX) 

a) 

Col 

46-50 

Vertical  incremental  distance  (INCY) 

(i) 

Third  Card 

Col 

1-5 

Print  mode  (IC0ST) 

a) 

Col 

6-10 

Number  of  blocked  nodes  (IBL0CK) 

(i) 

Col 

11-15 

Indicator  for  applying  interpolation  (INTER) 

(i) 

Col 

16-20 

Number  of  decimal  places  (NDEC) 

(D 

Col 

21-25 

Indicator  for  including  roughness  factor  (IR) 

(i) 

Col 

26-30 

Number  of  interpolated  data  points  +1  measured 

horizontally  (IGRIDX) 

a) 

Col 

31-35 

Number  of  interpolated  data  points  +1  measured 

vertically  (IGRIDY) 

a) 
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Fourth  Card 


Col 

1-5 

Below-sea-level  altitudes  (ZER0) 

(F) 

Col 

6-10 

c  N 

(F) 

Col 

11-15 

G 

(F) 

Col 

16-20 

H 

(F) 

Col 

21-25 

ALPHA 

}  Cost  function  parameters 

(F) 

Col 

26-30 

P 

1 

(F) 

Col 

31-35 

BETA 

| 

■(F) 

Grid  Altitude  Input  Data 

For  convenience,  the  model  is  programmed  to  accept  data  in  differ 
ent  formats  as  input  through  the  use  of  a  variable-format  statement. 
Fortran  provides  for  variable-format  statements  by  allowing  a  format 
specification  to  be  read  into  an  array  in  storage  and  using  the  infor¬ 
mation  in  the  array  as  the  format  specification  for  subsequent  input/ 
output.  Thus,  the  user  is  given  much  flexibility  in  the  form  of  the 
data  inputs. 


Fifth  Card 

Col  1-72  Variable  format  specification  of  the  input  data  (F0RMT)  (A) 
The  general  form  of  the  data  input  is  given  below: 


D 

D 


11 

21 


D 

D 


In 

2n 


ml 


m2 


D 

mn 


The  entries  D. .  in  the  table  are  the  data  values,  and  the  whole  array  of 
iJ 

these  numeric  values  is  called  the  input  data  matrix. 

The  general  data-card  input  is  keypunched  row-wise.  That  is,  all 
of  the  data  values  of  the  first  row  of  the  input  matrix  are  punched  in 
order  on  one  or  more  cards.  Then  starting  on  a  new  card,  the  second  row 
is  punched,  etc.  The  data  are  read  into  the  model  in  a  similar  fashion. 


Sixth  Cards 


!'8'  A-  P  4- 


Col  As  specified  by 

the  variable-  Input  data  values  (D.  .) 
format  ^ 


k 

Seventh  Card 


Col 

1-3 

X 

coordinate 

of 

first  node  to  be  blocked 

(i) 

Col 

4-6 

Y 

coordinate 

of 

first  node  to  be  blocked 

(i) 

Col 

7-9 

X 

coordinate 

of 

second  node  to  be  blocked 

(i) 

Col 

10-12 

Y 

coordinate 

of 

second  node  to  be  blocked 

(i) 

Col 

54-57 

X 

coordinate 

of 

tenth  node  to  be  blocked 

(i) 

Col 

58-60 

Y 

coordinate 

of 

tenth  node  to  be  blocked 

a) 

If  there  are  more  than  10  nodes  to  be  blocked,  continue  on  another  card 
of  the  same  form,  until  IBL0CK  nodes  have  been  specified. 


Dynamic  Program  Data 

After  creation  of  the  network  system,  the  model  will  predict 
either  optimal  routes  for  the  entire  network  or  some  subset  thereof,  as 
specified  by  the  user.  This  option  is  an  economy  consideration  in  that 
the  user  may  want  predicted  routes  for  different  areas  of  the  network. 
An  8-order  card,  defining  the  boundaries,  is  required  for  each  area. 


Eighth  Card 

Col  1-5 

X 

coordinate 

of 

lower  bound 

(ID0X) 

(I) 

Col 

6-10 

X 

coordinate 

of 

upper  bound 

(IT0X) 

(I) 

Col 

11-15 

Y 

coordinate 

of 

lower  bound 

(JD0Y) 

(I) 

Col 

16-20 

Y 

coordinate 

of 

upper  bound 

(JT0Y) 

(I) 

Col 

21-25 

Number  of  iterations  before 

test  (ISTART) 

(I) 

Col 

26-30 

X 

coordinate 

of 

terminal  node 

(I) 

Col 

31-35 

Y 

coordinate 

of 

terminal  node 

(I) 

k 

Include  only  if  IBL0CK  >  0. 
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i  As  many  sets  of  cards  1-8  as  there  are  cases  involving  different 
cost  calculations  may  be  included  during  a  computer  run,  with  each  set 
separated  by  a  blank  card*  The  end  of  an  input  deck  is  specified  by 
an  end-of-file  card. 


Appendix  C 

FLOW  CHARTS  OF  THE  PRINCIPAL  SUBROUTINES  TO  PREMIR 
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Main  Program 


Subroutine  BORDER 
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Subroutine  INTERP 


Calculate  squared  distance 
from  input  grid  points 
to  interpolated  data 
point  and  weight  each 
grid  data  point  accordingly 


Round  interpolated  value 
to  "NDEC"  places 


Alt  (I,  J)  =  rounded 
interpolated  value 


Alt  (IF,  JF)  =  -  999 
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Subroutine  COST 
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Subroutine  DYNAMI 
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Function  Subroutine  AMIN 


Yes 
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*? 


Appendix  D 

FORTRAN  IV  SOURCE  LISTING  OF  PREMIR 


on  o  ooo  ooo  ooo  non  ooo 


-70- 


COMMON  /ARRAYS/  COST ( 100, 100,8 ) , VALUE ( 1 00 , 1 00 , 2 ) , T ( 2  ) 

COMMON  /VALUES/  MTOl X , NT0 1 Y , NT0X , NT0Y , ND0X , ND0Y , I NCX , I NCY , I  OLD, 

1  I , J , I J J , NUMX , NUMY , I  COST, I  BLOCK, NDEC, IGRIOX , I  OR  I OY , Z FRO , C , G, H , 

2IST ART, ALPHA, P, BETA, IR, INTER 
DIMENSION  Z ALT ( 100, 100 ) ,  ALT (  100, 100 ) , FORMT ( 18 ) , A NAME ( 18 ) 
EQUIVALENCE  ( Z AL T ( 1 )  , COST ( 1 )  )  ,  < ALT (  1 ), VALUE (  1  >  ) 

INPOT  PROBLEM  CODE  NAME 

1  REA D( 5,503, EN0  =  60)  ( ANAME ( K )  ,  K  =  1 , 18 ) 

WR I TE ( 6, 506 )  ( ANAME (K) ,K=1 , 18 ) 

INPUT  MODEL  CONTROL  DATA 

I  F  R  R  0  R  =  2 

READ) 5 ,500, ERR  =  51 )  NUMX , NUMY , NUMX 2 , NUMY 2 , NI , MJ , ND0X,NI)0Y, I NCX , I NCY 
I  t  R  R  0  R  =  3 

READ ( 5 ,500, ERR  =  51 )  I  COST , I  BLOCK , I  NT ER , NO  EC , I R , I  OR  I DX , I  OR  I DY 
IF( IR.OT.O)  WRITE (6,600) 

INPUT  COST  DATA 


I  F  R  R  0  R  =  4 

READ ( 5 ,502 , ERR =51 )  Z ER 0 , C ,0 , H , AL PH A , P , BE T A 
I  E  R  R  0  R  =  5 

INPUT  OR  ID  ALTITUDE  DATA  FORMAT 

R  E  A  D  (  5 , 5 0 3  )  (  E UR  -i  T  (  I  )  ,  I  =  1  ,  1 8  ) 

WR ITE ( 6, 507 )  ( FURmT (  I  )  , I =1 ,  18  ) 

INPUT  GRID  ALTITUDE  DATA  INTO  ARRAY  ZaLT 

DO  10  J=NJ,uUMY 

R  E  A  0  (  5  ,  FOR  r-iT  ,  E R "'  =  5 0  )  (  Z  A  L  T  (  I  ,  J  )  ,  I  =N  I  ,  NUmX  ) 

10  CONTINUE 


C 


C 

C 

c 

c 

c 


THE  FOLLOWING  Ci  input F  0  VALUES  ARE  USED  IN  SUBSEQUENT  RuuTINFS 
EUR  CONTROLLING  DO-LOOPS,  FTC. 


NT  ox  =  I  GR  I  !>X  *  (  NUMX  2  —  l  )  +  NDOX 
N T I ) Y  =  I  G R  I  OY,  (  N Ui-i  Y 2 -  1  )  +  NDOY 

NT  0 1 X  =NT 0X+ I NCX 
,MTU1Y  =  MT0Y+  [NCY 

THE  FOLLOWING  Call  IS  FUR  SETTING  UP  BORDER  CONDITIONS 

call  border 


THE  FOLLOWING  Call  IS  EUR  INTERPOLATING  BETWEEN  INPUT  GRID  POINTS 


CALL  INTFRP 
C 

C  THE  FOLLOW  I  NO  C'.aLL  IS  FUR  COnPDT  I  ■'  0  THE  COSTS  OF  TRAVELING  TO 

C  ADJACENT  POINTS  and  ADJUSTING  FOR  BLOCKED  NODES. 

C 


°°°  oooon 
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CALL  COSTS 

c 

c  the  following  call  is  to  the  dynamic  recursion  routine 

c 

CALL  DYNAMI 
GO  TO  1 

SO  WRITE (6, 504)  I,J 
CALL  EXIT 

5  I  WRI  TE (6, SOS  )  I  ERROR 
SO  CALL  EXIT 
STOP 

SOU  FORMAT ( ISIS) 

SOI  FORMAT ( 1 6F  5 , 2  ) 

602  FORMAT (  16FS.0  ) 

503  FORMAT! 18A4) 

504  FORMAT !'  1 ',  1  ERROR  OCCURRED  WHILF  READING  INPUT  GRID  DATA.  COORDIM 

1  AT  E  VAt.UFS:  X  =  '  ,  I  4  ,  '  Y=',I4) 

505  FORMAT ('  1 ',  5X  ,' ERROR  OCCURRED  WHILE  READING  INPUT  CARD  NO.  ',13,'.' 
1  ) 

506  FORMAT  ( ' 1 ' , 2 X , ' PROBL EM  CODE :  '  ,  1 R A 4 ) 

507  FORMAT  {' 0 ',  2X  ,'  INPUT  FORMAT  SPECIFICATION':  1 , 2  X  ,  l  R  A  4  ) 

600  FORMAT  (' 0 ',  2X  ,' R  FACTOR  INCLUDED  ]m  COST  C  ALCUL  A 1  I  (INS  .  '  ) 

END 

SUHRDUT I  OF  BORDER 

COMMON  /ARRAYS/  COST ( 100 , 100, 6 )  , V  A  LOP (  100 ,100,2)  , T ( 2 ) 

COMMON  /VALUES/  NT  01 X  ,  N  (  0 1  Y  ,  NT  OX  ,  NTCY  ,  NOOX  ,  NOflY  ,  I  OCX  ,  I  M()Y  ,101.0, 

1  I  ,J,  I  J  J  ,  NUMX  ,  NlJNiY  ,  ICOST,  I  BLOCK  ,  NDFC  ,  IGRInX,  I  GR  I  DY  ,  Z.  FRO,C  ,  G,  H  , 
2ALPHA , P , HFT  A ,  I R 

DIMENSION  ZALTI 100,100) ,  A L T  (  I  00 , 1 Oo ) 

EDO  I  VALENCE  (  ZALT  (  1  )  ,COST.(  1  )  )  ,  (  A  L  T  (  1  )  ,  V  A  L  U  F  (  1  )  ) 

SUBROUTINE  ASSIGNS  SPECIAL  ALTITUDES  TO  BORDER  POINTS 
NT01X  IS  THE  UPPER  LIMIT  IN  THE  X  DIRECTION 
N 1 U 1 Y  IS  THE  UPPER  LIMIT  IN  THE  Y  DIRECTION 


N  I  =NDOX- 1 NC  X 
N J  =  NDO Y-  I  N C  Y 
DO  S  1 =  N I , N  T  0 1 X , 1 N  C  X 
ALT (  I  ,N.)  ) =-999 
ALT! I , NT01Y ) =-999 
S  CONTINUE 

DO  10  J=K'J,NT01Y,  INCY 
ALT (N I , J ) =-999 
ALT  (  NT01 X , J ) =-999 
It)  CONTINUE 
RETURN 
END 

SUBROUTINE  INTERP 

COMMON  /ARRAYS/  COST ( 1 00 , 1 00 , B  > , VALUE ( 1 00 , 1 00 ♦ 2 )  ,  T  (  2  > 

COMMON  /VALUES/  NT01 X , NTO 1 Y , NT OX , N TO Y , NDOX , NDOY , IMCX , I NCY, I  OLD, 
1 1  »  J  ,  I  J.J,  N1JMX  ,  NUMY  ,  ICOST,  I  BLOCK,  NO  EC,  I  GR  I  OX  ,  I  GR  I  D  Y  ,  Z  E  RO ,  C  ,  G ,  H  , 
2ISTART, ALPHA, P, BETA, IR, INTER 

DIMENSION  I F ( 10) , JF ( 10 ) 

DIMENSION  ZALT ( 100, 100) ,  ALT(100,100) 

EQUIVALENCE  ( ZALT ( 1 ) ,COST ( 1  )  )  ,  ( ALT ( 1 ), VALUE ( 1 ) ) 

INTERPOLATION  SCHEME 


on  20  I  =NOOX , NTOX , I NCX 


n  n  n  n  n  n 
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no  20  J=MD0Y  ,NT0Y t  INCY 
I  F  (  INTER.  GT.O)  0,0  TO  11 
A  =  (  I-2.*INCX)/IGRIDX 
R=  (  J-2 . *  I NC Y  )  /  I0,R  IOY 
I  A  =A 
1  h  =  H 
A=A-  I  A 
H=B— I R 

I  F  (  A  +B  .  F  0 . 0  )  GO  TO  12 
GX=0 
0  X  =  0 

NST1=IA+1 
NST2=IB+1 
NST3=NST  1  +  3 
NST4=NST2+3 

I F ( NST3. GT. NUMX )  NST3=NUMX 
I  F  (  NST4.GT  .NIIMY  )  NS1  4  =  K'IIHY 
C 

C  CALCULATE  SO.  HIST.  FROi‘1  INPUT  grid  POINTS  TO  INTERPOLATION 

C  POINTS  ARM)  WEIGHT  EACH  GRID  nAl  a  POUT  aCCURUINGIY 

C 

00  10  I  X  =  N S T  1  ,  N  S  T  3 
OU  10  I Y=NST2 , HST4 

0 1  S  T  2  =  (  (  I  -  2  .  *  I  R  C  X  )-IGKIOX*(  IX-2.  )  )**?.+(  (  J  -2  .  -  1  NC  Y  )  -  I  GK  I  DY*  (  IY-2.  ) 
A  )  *  *  2  . 

I)X=  0X+  1./0IST2 
GXA=  GX+  Z  A  L  T  (  IX.  I  Y  )  / 0 1  S T 2 

10  CONTINUE 

USE  *10'  TO  ROUND  OFF  INTERPOLATED  VALUE  TO  'NOEC  0 F C INAL  POINTS. 

I0=(  GX/DX  +  (  3  .  )  *(  10.  )  =**  (  -  (  MDEC  +  1  )  )  )*10**N0PC 

ALT  (  1  ,  J  )  =  I  (0  /  (  1  v  .  ( ) )  4  M  D  F  C 
GO  TO  20 

11  I  A  =  1-2 
I  P  =  J  -  2 

INTERPOLATION  IS  AT  AN  INPUT  DATA  POINT  UR  INTERPOLATION  IS  NOT 
BFING  USED. 

12  ALT(I,J)=  Z ALT ( I A  +  2  » I B+2 ) 

20  CONTINUE 

I  F (  I  BLOCK . L  E . 0 )  GO  TO  30 
I  FL AG  =  0 
B  1  =  I  BLOCK/ 10. 

I  F  =  H  1 
B2  =B 1-  I  B 
NUM=B1 

I  F ( B 2 . G T . 0 )  NUM  =  MUM+1 
22  I  FLAG- I FLAG+1 

R FAD (5,500)  (  I  F ( K )  ,  J F ( K ) , K  =  1 , 1 0 ) 

K= 1 0*B2 

IF(  (Bl.GE.l) . A  N  0 . (  IFLAG.LE.B1)  )  K  =  10 
OD  2  3  I K  =  1 , K 
ALT (  IF (  IK  )  ,  JF (  IK )  ) =  -999 
25  CONTINUE 

I  F (  IFLAG.LT. NUN)  GO  TO  22 
30  RETURN 
500  FORMAT (2013) 


o  o  o  o  o  o  o  o  o  o  o  o  o  O  o  o 
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ENO 

SUBROUTINE  COSTS 

COMMON  /ARRAYS/  COST  (  100, 1 00  ,  ft  )  ,  VALUE  (  100, 100,2  )  ,T(  2  ) 

CO  I'm  MON  /VALUES/  NT01X  ,NT01  Y  ,  NTOX  .  M  TOY  ,  NOOX  ,  MOOY  ,  INCX  ,  I  MCY  ,  I  OLD, 

1  I , J , I J J , NIWX,  NUMY , I  COST .  I  BLOCK, NOEC , I  OR  I  OX , I  OR  I OY , Z  E RO , C , G , H , 

21  START  ,ALPHA,P , BETA, IR , INTER 
DIMENSION  7  ALT (  100, 100 )  ,  ALT ( 100, 100) 

EOU I  VALENCE  ( Z A L T (  1  )  ,COST II)),  ( ALT (  1 )  ,  VALUE  ID) 

DO  20  I=NOOX,NTOX, INCX 
00  20  J  =NDOY , N TOY , I MCY 

EIGHT  WAYS  OUT  OF  EACH  I , J  POINT 

00  2  0  I JOUT  =  1 , R 
I  T 0= I +GOTOX ( IJOUT) 

JTO- J+GOTOY ( IJOUT ) 

I  F (  ALT( ITO,  JTO) .LE. ZERO)  GO  TO  10 
A  =  1  JOUT / 2 . 

I  A  =  A 
A  =  A  -  I  A 

SLOPE= ( H/G ) *(  ALT ( ITO, JTO l-ALT ( I , J )  ) 

IE(A.EO.O)  GO  TO  5 

COST  CALCULATION  FOR  TRAVELING  TO  A  MON-U I  AGON  A I.  POINT 
COST ( I , J , I JOUT ) =  <  1.  +  C*( ARS(aLPHA*SLorE) )*#P  ) 

I  F  (  IR  .GT.  0  )  COST  (  I  ,  J  ,  I  JOUT  )  =CUST  (  I  ,  J.  I  JOUT.)  +  (  0  E  T  A  *  (  A  L  T  I  I  T  0+  1  ,  J  T  O 

1  )  -  ALT  (  I  T  0-1  *  JTO  )  )  /  (  2*G  )  )  **2  +  I  B  ETA*  (ALT  (  IT  O  .  J  T 11+  1  )  -  AL  T  (  I  TO,  JTO- 1  ) 

2  )  /  (  2  *  G  )  )  *  *  2 
GO  TO  20 

COST  CALCULATION  FOR  TRAVELING  TO  a  IUanOMAL  POINT 

5  COST  (  I  ,J,  IJOUT  !=  (  1.4  +  C*(  ARS(  .  7*  A  I.  Ph  A*  SLOP  E  )  )**P) 

I F (  IR.GT.O)  COST (  I  .  J,I JOUT ) =COST (  I , J, I JOUT )  +  (  B  E  T  A  *  (  A  |_  T  (  I T  0+ 1 , J  T  0 

1 )  -ALT  I  I TO-1 , JTO  )  )  / ( 2*G )  ) **2  +  ( BETA* ( ALT ( T  T  0 , J T 0+ 1 ) - A L T ( I  TO, JTO- 1 ) 

2 )  / ( ?*G)  ) *  *  2 
GO  TO  20 

ATTATCH  A  HIGH  COST  TO  POINTS  MTH  ALTITUDES  LESS  THAN  OR  EOU  AL  TO 
THE  VARIABLE  'ZERO'. 


10  COST ( I , J, IJOUT ) =10. **6 
20  CONTINUE 

IF  '  I.COST  '  IS  GREATER  THAN  ZERO,  PRINT  OUT  COST  CALCULATIONS 

I  E ( I  COST . L  F .0 )  GO  TO  40 
CALL  PAGE (LIME) 

DO  30  I=NDOX , NTOX , INCX 
DO  30  J  =  N DO Y , N T 0 Y  ,  1 N C Y 
L INE=L INE+1 

IF(LINE.GT.55>  CALL  PAGF(LIME) 

WR ITE ( 6,501 )  I , J , ALT ( I , J ) . ( COST ( I , J , I J ) , I J=1 , « > 

30  CONTINIIF 
40  RETURN 

501  FORMAT ( IX, 14, 15 , E9. 1 ,9X,RG1 1.6  ) 

END 

SUBROUTINE  OYNAMI 


ono  ooo  non  non 
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COMMON  /ARRAYS/  COST  (  100, 100,8  )  ,  VALUE  (  100, 100, 2  )  ,T(  2) 

COMMON  /VALUES/  NT  0 1 X  ,  NT  01 Y  ,  NTOX  ,  NTflY  ,  MUOX  ,  NDOY  ,  INCX,  INCY,  I  OLD, 

1  1  ,  J  ,  I  J  J  ,  NIJMX  ,  NUMY  ,  I  COST  ,  1  BLOCK  ,  NO  EC  ,  I  OK  I  OX  ,  IGR  1  DY  ,  ZERO  ,  C  » 0  ,  H  , 

2  I  START, ALPHA, P, BETA , IK,  INTER 
DIMENSION  POLICY ( 100, 100) 

INTEGER  T 

INTEGER*?  POLICY 

3  R  E  A  0 ( 5,515, END=60 )  I  OOX , I TOX , J OOY , J TO Y ,  I S  TAR  T , T (  1  )  ,T(2) 

IF( IDOX.LE.O)  RFTURN 
IE( ISTART. LE.O)  ISTART=1 
00  5  I T  =  1 , 2 
do  5  i=innx,iTnx 
DO  5  J= JOOY  ,  JTOY 
VALUE (  I  ,J, I T ) =0 

5  CONTINUE 

1  K'  IS  THE  ITERATION  NUmBFR 
K  =  0 

6  K=K+ 1 

CALL  F I  X (  I  NEW,  I <  L  H , K ) 

DYNAMIC  MINIMIZATION  SCHEME 

DO  10  I =IOOX, I TOX , INCX 
00  10  J  =  JOOY , JTOY , INCY 

IF( (  I .EO.TI 1 )  )  .AMO. ( J.E0.T(2 ) )  )  GO  TO  10 
VALUE ( I  , J , I  NEW ) =AMI N ( I OOX , I TOX,  JOOY , JTOY ) 

POLICY!  I  ,J)=  I J J 
10  CONTINUE 

DO  'ISTART'  ITERATIONS  BEFORE  CHECKING  LIST  FOR  TERMINATION. 

I  F (K. LE . ISTART )  GO  TO  6 
DO  15  I  =  I  00 X  ,  I T f  =  X , INCX 
DO  15  J=JDOY, JTOY, INCY 

IF!  VALUE (  I  ,  J  ,  I  NEW ) .ME . VALUE (  I , J , I  OLD)  )  GO  TO  6 
15  CONTINUE 

RECURSION  TERMINATED;  OUTPUT  OF  RESULTS  FOLLOWS. 

WRITE (6,500)  K  ,  T (  1  ) , T ( 2 ) 

L  I  N  E  =  0 

DO  20  I  =  IDOX,  ITOX, INCX 
DO  20  J  =  JOOY,  JTOY, INCY 

IF!  (  I.EO.T(l)  )  . A  NO . ( J . F  o . T ( 2 )  )  )  GO  TO  19 

L I NE=L INE+1 

IF(LINE.GT.52)  WRITE (6, 500)  K,T(1),T(2) 

IF(LINE.GT.52)  LINE=1 
I JOUT=POL ICY ( I , J ) 

I P= I +GOTOX ( IJOUT) 

J  P  =  J  +  GOTOY (  IJOUT) 

WRITE(6, 510)  I,J,  IP,JP,  VALUE ( I , J, INEW) 

GO  TO  20 

19  WRITE! 6,520)  I  ,J 

20  CONTINUE 
GO  TO  3 

60  CALL  EXIT 

500  FORMAT!  ' 1 ' ,4X,  'NUMBER  OF  ITERATIONS  REDO  I  RED  TO  REACH  SOLUTION:', 


onno  o  o  o  ooooo 
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1 1  6 , // t 5X, 1  FROM  POINT  GO  TO  POINT  MIN  COST  TO  TFRMINAL ( X= 1 , 

2  13,  ' ,  Y=', 13, ')'»/, 

3  7X,'X  Y ' , 9X  ,  '  X  Y ',/, IX, 1301  «*•■),/ ) 

510  FORMAT (2X,2I6,5X,2I5»  5X,F20.2> 

515  FORMAT ( 1615) 

520  FORMAT ( 2X ,2 16 , 5X ,' STOP — TERMINAL  POINT.') 

RFTURN 

END 

SUBROUTINE  F I  X ( INEW,IOLD,K) 

SUBROUTINE  GIVES  ALTERNATING  VALUES  OF  '1'  OR  TO  INEW  AND  IOLD 

DEPENDING  ON  WHETHER  K  IS  EVEN  OR  ODD:  K  EVEN  GIVFS  INFW  VAIUF 

OF  '2'. 

5  A  =  K  /  2 . 0 
I  A  =A 
A  =  A  -  I  A 

IF  A  IS  ZERO  K  IS  EVEN 

IF(A.FO.O)  GO  TO  10 
I  N  E  W  =  1 
I OL  D  =  2 
RETURN 
10  I N  E  W  =  2 
I  0  L  0  =  1 
RETURN 
END 

FUNCT  I  UN  AM  I  N  (  I  DUX  ,  I  T  OX  ,  J  UilY  ,  J  TOY  ) 

COMMON  /ARRAYS/  COST (  100 , 1 OO , 8 ) , VALUE (  100,100,2)  ,  T (2) 

C C i M MOW  /VALUES/  WTO  1  X  ,  NTO 1  Y  ,  NTOX  ,  NTOY  ,  NDOX  ,  N D 0 Y  ,  I  w C X  ,  I  N C  Y  ,  I  U L D  , 
11*0, I J  J , NUMX , NUM Y , I  COST , I  BLOCK , NDFC ,  IGRIDX, I GR I DY , Z FRO , C , G , H , 

21  START, ALPHA, P, BETA. IR,  INTER 
INTEGER  T 

SUBROUTINE  FINDS  A  RELATIVE  m I N I  HUM  VALVE  OUT  OF  EACH  POINT  AND 
ASSIGNS  THE  VALUE  T'~  AMIN.  THE  PATH  POINT  IS  ASSIGNED  TO  IJJ. 

A= 10**25 
DO  5  I  T  =  1 , 8 
I P  = I +  GOTOX (  IT) 

J  P  =  J+GOT OY (  IT ) 

I F (  ( I P. L T . I  DUX ) . OR .  (  I  R . GT . I TOX )  )  GO  TO  8 
IF ( ( JP.LT. JDOY ) .OR . ( JP.GT. JTOY) )  GO  TO  5 
I F (  VALUE (  IP, JP , I OL  1  +  COST ( I , J, I T )  . L T . A  )  GO  TO  4 

GO  TO  8 

4  A  =  VALUE (  IP, JP, IOLD)  +  COST (1,0,  IT) 

I JJ=IT 

5  CONTINUE 
AN  I N  =  A 

R  r  T  i  )R  N 
END 

SUBROUT  I ME  PAGE ( L I  WE ) 

COMMON  /ARRAYS/  COST! 100, 100, « ), VALUE ( 100 , 100, 2 ) , T ( 2 ) 

COMMON  /VALUES/  NT 0 IX, NTO 1Y, NTOX, NTOY, NDOX , NDOY , INCX, INCY, IOLD, 

II  ,J,  I  J  J  ,NIJMX  ,  Nl  IN  Y  ,  IC'iST  ,  I  BLOCK  ,  NDFC  ,  IGRIDX,  I  GR  I  DY  ,  Z  ERG  ,  C  ,  G  ,  H  , 

2  I  START, ALPHA , P , BETA . I R , I NTER 

SUBROUTINE  PRODUCES  a  PAGE  HEADING  WHEN  OUT  PUTTING  COSTS  RESULTS. 
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C 

C 

C 


c 

c 

c 


c 

c 

c 


/* 

// 


INTEGER  T 
LINE  =  0 

WRITE (6,500)  INCX, INCX, I MCY , I NC Y , I NCX ,  I  NCY ,  INCX, I NCX , I NCY, I NCY, 
1 INCX, I NCY 
RETURN 


500  FORMAT  {  *  1  *  ,4‘OX  ,  ***COSTS  OF  TRAVELING  FROM 

1  *  1  * / / ♦ 3 0 X  t •POINT  1  POINT  2  POINT  3 

2  P  0  I N  T  6  POINT  7  POINT 

3  4X  t  *  X  y  ALTITUDE  COSTS:  X+',I2,»,  V  ; 

4  ,  Y  +  1  , 1  2  ,  '  :  x-'  ,I2t  1  tY+«  ,12,  1  ;  X-',I2,*,  Y 
5X  ,  Y-  *  ,  12  ,  1  ?  X  +  ■  ,  12  ,  '  ,  Y-  1  ,  I  2,  *  .  1  ,  /  ,  IX  ,  130  ( 

END 


X , Y  TO  ADJOINING  POINTS* 
POINT  4  POINT  5 

X+« , 12, • ,Y+» , 12, 1 :  X 

:  X-' , 12, 1 ,Y-* , 12, 1 ; 

1  *  *  )  ) 


FUNCTION  GOTO ( DUMMY  ) 


COMMON  /VALUES/  NT01X  ,.NT  01  Y  ,NTOX  ,  N  TOY  ,  NO  OX  ,  MDOY  ,  I  NCX  ,  I  NCY  ,  I  OLD, 

1  I  ,  J  ,  I  J  J  ,  N  U  MX,  NUMY  ,  I  COST  ,  I  BLOCK  ,  Mf)FC  ,  I  GR  I  r»X  ,  I  G  R  I  D  Y  ,  Z  E  R  0 ,  C  ,  G ,  H  , 

2  I  START  t  A  L  PH A , P  ,  BETA, IR  , INTER 


SUBROUTINE  IS  USED  FOR  FINDING  DIRECTION  OF  THE  OUT  POINTS. 


G  O  T  0  =  0 
RETURN 


OUT  POINT  IN  THE  X  DIRECTION 
ENTRY  GOTOX ( I JOUT  ) 

G  n  t  n  x = - 1  n  c  x 

I  F  (  {  I  JOUT  . LT . 3 ) . OR . ( I  JOUT  .  FO.  8  )  )  GOTOX  =  I  CX 
IFI  (  I  JOI JT  .  FD  .  3  )  .  OR  .  (  I  JOUT  .  ED  .  7  )  )  GOTOX  =  0 

R  tr  T I  j  R  M 


OUT  POINT  IN  THE  Y  1)1  RUCTION 

1 0 TRY  GOTOY(IJOOT) 

G  0  T  n  Y  =  -  I  i\)C  Y 

I  F  (  (  I  JOMT  .  PC).  1  )  ,0R  .  (  I  JOUT  .  FO.  5  )  )  GOTOY  =  0 
I  F  (  (  I  vliil  'T  .  GT  .1  )  .  A  Kill,  (  I  JOUT  .  L  1  .  S>  )  )  GOTOY =  !  -MCY 
R  E  T  0  R  w 
E  M  0 
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